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1 QCD Sum Rules: Twenty Years After 



I will discuss a method of treating the nonperturbative dynamics of QCD which 
was created almost twenty years ago in an attempt to understand a variety of 
properties and behavior patterns in the hadronic family in terms of several basic 
parameters of the vacuum state. The method goes under the name QCD sum rules 
- rather awkward, for many reasons. First and foremost, it does not emphasize 
the essence of the method. Second, in Quantum Chromodynamics there exist many 
other sum rules, having nothing to do with those suggested in Ref. Finally, 
some authors add further confusion by using ad hoc names, e.g. the Laplace sum 
rules, spectral sum rules, and so on, which are even foggier and are not generally 
accepted. 

It would be more accurate to say "the method of expansion of the correlation 
functions in the vacuum condensates with the subsequent matching via the dis- 
persion relations". This is evidently far too long a string to put into circulation. 
Therefore, for clarity I will refer to the Shifman-Vainshtein-Zakharov (SVZ) sum 
rules. Sometimes, I will resort to abbreviations such as "the condensate expansion" . 

Twenty years ago, next to nothing was known about nonperturbative aspects of 
QCD. The condensate expansion was the first quantitative approach which proved 
to be successful in dozens of problems. Since then, many things changed. Various 
new ideas and models were suggested concerning the peculiar infrared behavior 
in Quantum Chromodynamics. Lattice QCD grew into a powerful computational 
scheme which promises, with time, to produce the most accurate results, if not for 
the whole set of the hadronic parameters, at least, for a significant part. 

It seems timely to survey the ideas and technology constituting the core of the 
SVZ sum rules from the modern perspective, when the method became just one 
among several theoretical components in a modern highly competitive environment. 
An exhaustive review of a wealth of "classical", old elements of the method and 
applications was given in Ref. 0|. There is hardly any need in an abbreviated 
version of such a report. New applications which were worked out in the last decade 
or so definitely do deserve a detailed discussion. As far as I know, no comprehensive 
coverage of the topic exists in the literature. Unfortunately, in these lectures I will 
not be able to provide such a coverage, which thus remains a task for the future 
Q. Instead, I will focus on those qualitative aspects where understanding became 
deeper. This is the first goal. Secondly, selected new applications will be considered 
to the extent that they illustrate the theoretical ideas of the last decade. And last 
but not least, I will try to outline an ecological niche which belongs to the SVZ 
method today. As a matter of fact, over the years, slow but steady advances were 
taking place in our knowledge of the hadronic world. Some old and largely forgotten 

"'^Work on systematically reviewing a variety of developments that took place since the mid- 
1980's and numerous new applications is under way (a private communication from B.L. loffe). 
A survey devoted to the relation between the sum rule and lattice results is being written by A. 
Khodjamirian. 
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predictions of the SVZ sum rules were confirmed recently by other investigations 
based on totally different principles. These predictions are extremely nontrivial. 
For instance, about 15 years ago, it was discovered that not all hadrons are alike; 
there are remarkable distinctions between them, especially in the glueball sector. 
The fact that not all hadrons are alike is now becoming more and more evident 
from the lattice results as well. Other examples of this type are known too. By 
confronting them with alternative sources of information, such as lattices, we get 
a much fuller picture of the QCD vacuum. This process might be very beneficial 
to both sides. Unfortunately, at present the lattice and analytic QCD communities 
are largely disconnected, and rarely talk to each other. My task is to show that 
borrowing from each other makes everybody richer! 



2 QCD Vacuum and Basics of the SVZ Method 
2.1 General ideas 

The color dynamics described by QCD is very peculiar. If we have two probe color 
charges, their interaction approaches the Coulomb law at short distances, with a 
weak coupling constant. The Coulomb interaction is due to the one gluon exchange 
(Fig. 1). At larger distances, the gluon starts branching (Fig. 2), which leads to 
a remarkable phenomenon known as antiscreening, or asymptotic freedom Q]. In 
normal theories, like QED, the virtual cloud screens the bare charge making the 
charge seen at larger distances smaller than the bare one. This situation is perfectly 
transparent intuitively. 

In QCD, instead of screening, the branching processes result in a totally counter- 
intuitive behavior - the antiscreening. Those of you who would like to know the 
physical origin of antiscreening are referred to the very pedagogical review |^, Sect. 
1.3. If the distances are not too large one can apply perturbation theory to quanti- 
tatively describe the gluon branchings, and derive the famous formula of asymptotic 
freedom 

as = Const/ In — (1) 
r 

where tq is a dynamically generated scale parameter of QCD, rg ~ 1 fm. At small 
separations the effective coupling constant dies off logarithmically. 

Usually one considers the running coupling constant in the momentum space. 
Then in the leading (one- loop) approximation corresponding to Eq. (|l]), 

where is the number of colors and Nf is the number of flavors and b is the first 
coefficient in the Gell-Mann-Low function. In the low-energy domain we will be 
interested in (below the threshold of the charm production) the number of active 
flavors is Nf = 3, and, hence, 6 = 9. 
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Figure 1: Heavy color charges (solid lines) interacting through the one-gluon ex- 
change (dashed line). 



In perturbation theory the effective coupling is calculated order by order. The 
result is known up to three loops, and for Nc = Nf = 3 can be conveniently written 
as H 



/ 2x 1 Inln^ (0.79)2 

a = -0.79 , f + ^ , i 

' In^ In^^ In^i^^ 



A2 A2 A2 



2 2 

(lnln^)2-lnln^ + 0.415 

^ IS? IS? 



+0 



1 



(3) 



where we define 



(4) 

4 TT 

and A is the scale parameter of QCD introduced in a standard way 0. In the third 
and higher loops the law of running of becomes scheme-dependent. The third 
term in Eq. (^ refers to the so called modified minimal subtraction (MS) scheme. 
More exactly, since Eq. (^) describes running with three flavors, the parameter A 

is actually A^^^^. Below we will deal exclusively with Ai^Tr; therefore, not to make 

MS MS 
the notation too clumsy, we will suppress the sub (super) scripts. 

Let us leave a while the issue of the effective coupling constant, with the intention 
of returning to it later. The only lesson one should remember at this stage is that 
the effects caused by perturbative gluon exchanges are logarithmic. 

Equations (|I|) or (^ imply that the effective color interaction becomes stronger 
as the separation between the probe color charges increases. Being remarkable by 
itself, this phenomenon carries the seeds of another, even more remarkable property 
of QCD. When the distance becomes larger than some number times A~^, and 
exceeds a critical one, the branchings of gluons become so intensive (Fig. 3) that it 
makes no sense to speak about individual gluons. 

Rather, we should phrase our consideration in terms of the chromoelectric and 
chromomagnetic fields. In "normal" theories, like QED, the field induced by the 
probe charges separated by a large distance, is dispersed all over space. In QCD 
it is conjectured that a specific organization of the QCD vacuum makes such a 
dispersed configuration energetically inexpedient |^. Rather, the chromoelectric 
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Figure 2: Gluon branching responsible for the logarithmic running of the effective 
gauge coupling constant in perturbation theory (shown is a sample one-loop graph). 




Figure 3: When the separation between the probe color charges exceeds critical the 
nonlinearity becomes so strong that it makes no sense to speak about individual 
gluons. The color fields between the probe charges form a flux tube with transverse 
dimensions ~ A~^. 
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field between the probe charges squeezes itself into ge-like configuration. The 

situation is reminiscent of the Meissner effect in superconductivity. As is well-known, 
the superconducting media do not tolerate the magnetic field. If one imposes, as 
an external boundary condition, a certain flux of the magnetic field through such 
a medium, the magnetic field will be squeezed into a thin tube carrying all the 
magnetic flux; the superconducting phase is destroyed inside the tube, and there is 
no magnetic field outside. 

Superconductivity is caused by condensation of the Cooper pairs - pairs of elec- 
tric charges. A phenomenologically acceptable picture of the infrared dynamics in 
QCD requires chromoelectric flux tubes, not magnetic ones. If something like con- 
densation of "chromomagnetic monopoles" took place in the QCD vacuum then this 
would naturally explain nonproliferation of the chromoelectric field from the exter- 
nal probe color charges in the entire space - a dual Meissner effect would force the 
field to form flux tubes in the QCD vacuum, in this way ensuring color confinement 

i- 

Nobody ever succeeded in proving that the dual Meissner effect does indeed take 
place in QCD. There are reasons to believe, however, that it does. First, some 
evidence came from the numerical study in the lattice QCD (see e.g. p). Quite 
recently, a breakthrough was achieved in understanding of the infrared dynamics 
in a theory which might be considered a relative of QCD, the so called N = 2 
super symmetric (SUSY) Yang-Mills theory, in which it was analytically shown that 
the monopoles do condense in the strong coupling regime, thus providing a basis for 
the dual Meissner mechanism. 

Even if, as a result of a breakthrough sometimes in the future, it becomes clear 
how a similar mechanism might develop in QCD, we have a long way to go from a 
qualitative picture of the phenomenon to a quantitative approach allowing one to 
exactly calculate all the variety of the hadronic properties from the first principles. 
It may be that such string-based exact approach will never emerge. While the issue 
of the exact solution is still under investigation, let us see what analytic QCD can 
offer as an approximate solution. 

The basic idea lying behind the SVZ method is quite transparent. The quarks 
comprising the low-lying hadronic states, e.g. classical mesons or baryons, are not 
that far from each other, on average. The distance between them is of order A~^. 
Under the circumstances, the string-like chromoelectric flux tubes, connecting well- 
separated probe color charges, hardly have a chance to be developed. Moreover, 
the valence quark pair injected in the vacuum, in a sense, perturbs it only slightly. 
Then we do not need the full machinery of the QCD strings, whatever it might 
mean, to approximately describe the properties of the low-lying states. Their basic 
parameters depend on how the valence quarks of which they are built interact with 
typical vacuum field fluctuations. 

It is established that the QCD vacuum is characterized by various condensates 
m. Half-dozen of them are known: the gluon condensate Gf^^^, the quark condensate 
qq, the mixed condensate qaGq, and so on. The task is to determine the regularities 
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and parameters of the classical mesons and baryons from a few simple condensates. 

Of course, without invoking the entire infinite set of condensates one can hope 
to capture only gross features of the vacuum medium, at best. Correspondingly, any 
calculation of the hadronic parameters of this type is admittedly approximate. Since 
hadronic physics is deprived of small parameters, in the vast majority of problems 
high accuracy is by far not the most desired requirement to calculation, however. 
Does it make any difference if, say, the proton magnetic moment is predicted theo- 
retically to two or to three digits? I do not think so. Rather, it is high reliability of 
predictions in a wide range of problems where the answer is not known a priori and 
theoretical control over qualitative aspects that are the primary goals of the theory 
of hadrons. Below I will try to demonstrate that the SVZ method is both a reliable 
and controllable approach which has enjoyed enormous success in dozens and dozens 
of instances. There are a few cases where it fails 0, but we do understand why. 
And the very fact of the failure of the standard strategy teaches us a lot. This hap- 
pens for specific, "nonclassical" hadrons, with a very strong coupling to the vacuum 
fluctuations, which are very different from say the p meson or nucleon. So, we can 
reverse the argument and convert the failure into success by saying that the method 
predicts that not all hadrons are alike. 

2.2 Getting started /Playing with toy models 

The basic microscopic degrees of freedom of Quantum Chromodynamics are the 
quarks and gluons. Their interaction is described by the Lagrangian 

^ = -W,uG% + T.1^Pq (5) 
/ 

where / is the flavor index, and the color index of the quark fields q is suppressed. 
For simplicity we will assume that the quark mass terms vanish. In other words, 
we will work in the chiral limit. In this limit the pion is massless, = 0. This is 
known to be quite a good approximation to the real world. 

Neither quarks nor gluons are asymptotic states. Experimentally observed are 
hadrons - color-singlet bound states. In order to study the properties of the classical 
hadrons it is convenient to start from the empty space - the vacuum - inject there 
a quark- ant iquark pair, and then follow the evolution of the valence quarks injected 
in the vacuum medium. The injection is achieved by external currents. The most 
popular are the vector and axial currents. Their popularity is due to the fact that 
they actually exist in nature: virtual photons and W bosons couple to the vector 
and axial quark currents. Therefore, they are experimentally accessible in the e^e~ 
annihilation into hadrons or hadronic r decays. 

Thus, the objects we will work with are the correlation functions of the quark 
currents. More concretely, let us consider the vector current with the isotopic spin 
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1 = 1, 

= , or = wy^d. (6) 

The first current shows up in the e'^e~ annihilation while the second is relevant to 
the r decays. The two-point function U^^ is defined as 

n^. = ^/ e^''Vx(0|T{J^(a;) Ji(0)}|0) . (7) 

where q is the total momentum of the quark-antiquark pair injected in the vacuum. 
Due to the current conservation U^iy is transversal and, hence, 

^^.u = [q^^qu - q^g^^uMq^) ■ (8) 

For historical reasons n(g^) is often called the polarization operator, we will use 
this nomenclature in what follows. It is easy to count that the function Il{q^) is 
dimensionless. The imaginary part of n(g^) at positive values of (i.e. above the 
physical threshold of the hadron production) is called the spectral density, 

127r 

p{s) = —lmU{s), s^q\ (9) 



Up to normalization it coincides with the cross section of e^e~ annihilation into 
hadrons (measured in the units a{e~^e~ — > fi'^p")) or the r decay distribution func- 
tion. The numerical factor in the definition of p{s) in Eq. (^ is introduced for 
convenience, as will become apparent shortly. In the real world Nc = 3 but we will 
keep this factor explicit for a while, to keep track of the Nc dependence. 

The spectral density carries full information about the spectrum and widths of 
hadrons with given quantum numbers. Every QCD practitioner dreams of the exact 
calculation of the spectral density. Later on we will see how the SVZ sum rules 
constrain the parameters of the lowest-lying state, the p meson in the case at hand. 
But first, prior to submerging into technical aspects of the SVZ approach, let us 
make an educated guess of how the spectral density might look, in gross features, 
to get an idea of what is expected for n(g^). 

To begin with, consider the limit Nc —>■ oo. As is well-known ||10|, in this limit 



all hadrons are infinitely narrow, since all decay widths are suppressed by powers of 
l/Nc- Correspondingly, p(s) is a sum of delta functions. 

Moreover, there are good reasons to believe that these delta functions must 
be approximately equidistant, at least asymptotically, for highly excited states. A 
string-like picture of color confinement naturally leads to (approximately) linear 
Regge trajectories, and, hence, oc n at large n where m„ is the mass of the n-th 
p-meson excitation. In the world of purely linear Regge trajectories there are an 
infinite number of daughter trajectories associated with each Regge trajectory. The 
daughter trajectories are parallel to the parent trajectory and are shifted by integers. 
Thus, in the old Veneziano model (a review is given e.g. in |jll|) = + n/a' 
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Figure 4: The spectral density in the toy model. For clarity I gave a tiny width to 
the S functions in Eq. (|TU|). 



[a' is the slope of the Regge trajectory). In some of the later versions, which are 
more "QCD-friendly", = + 2n/a' , see e.g. [|12|. Since for now the focus is on 



a qualitative picture, the distinction between these two scenarios is not essential for 
our illustrative purposes. For definiteness, let us accept that the distance between 
two consecutive excitations contributing to p{s) is 2/a' ~ 2 GeV^. 

Assembling all these elements together one obtains a sketch of the spectral den- 
sity (Fig. 4), 

oo 

p(s) = 3 ^(5(s-l-3n). (10) 

n=0 

Here is set equal to unity; all dimensional quantities are measured in these units, 
for instance, 2 /a' ~ 3. The couplings of all mesons to the current are chosen to 
be equal. This choice is represented by the overall numerical factor 3 in the sum 
(PUp. The explaination for this will become clear momentarily. 

Now, given the imaginary part ([10|), it is not difficult to reconstruct the polar- 
ization operator H itself, 

n(g') = -^V^(^) + Const., (11) 

Q = -g , z = — - — , 



9 



where ip is the logarithmic derivative of the T function, 



k=0 



k + 1 z + k 



The constant on the right-hand side of Eq. ( [TTl) is irrelevant, since it can be always 
eliminated by an appropriate subtraction. It will be omitted hereafter. 

Although our desired target is the spectral density on the physical cut, i.e. at 
positive values of g^, all theoretical calculations in Quantum Chromodynamics are 
carried out off the physical cut, say, at negative values of (positive Q^). The 
reason is obvious: the QCD Lagrangian (^ is formulated in terms of quarks and 
gluons, not hadrons. In the absence of the final solution of QCD we can deal only 
with the quark and gluon fields. Working in the Euclidean domain, off the physical 
cuts, we can calculate in terms of quarks and gluons. This is a common feature of 
the SVZ sum rules and lattice strategies. 

At positive values of an asymptotic representation exists for the ip function, 

1 °° R 

(12) 

where i?2n stand for the Bernoulli numbers, 

i?2n = (-ir^S^C(2n); (13) 



(27r, 



here ( is the Riemann function. (In some textbooks (— l)"+^i?2n is called the n-th 
Bernoulli number and is denoted by Bn-) 

Equations ([TT|) and ([T2|) at large positive imply that in the leading approxi- 
mation 

n(Q')--T^inQ^ (14) 



This logarithmic formula for the polarization operator exactly matches what we 
expect from perturbation theory. Indeed, at large Q^, in the deep Euclidean domain, 
the points of injection and annihilation of the quark pair in the vacuum are separated 
by a small space-time interval. The quarks have no time to interact with the vacuum 
medium. They propagate as free objects (Fig. 5). For free quarks, obviously, 

= ^ J e"^^d'xTT{^^So{x,Oh,So{0,x)} , (15) 

where So{x,y) is the free-quark Green function describing propagation from the 
point y to the point x, 

1 ^ 

Soix^y) = ^-(^^ ^ = x-y. (16) 
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Figure 5: The two-point function of the vector currents in the free-quark approxi- 
mation. The external current injecting the quark- ant iquark pair in the vacuum (and 
then annihilating it) is denoted by wavy lines. 



It is trivial to obtain 



Tr {^,So{x, 0)7.^o(0, ^)} = ^'^'''^s "^^^"^ ■ (1^) 



Note that the expression on the right-hand side is automatically transversal. Sub- 
stituting Eq. {\r7\j in Eq. ([15D and doing the Fourier transformation, we arrive at 
Eq. dg). 

Thus, the leading asymptotic behavior of the polarization operator stemming 
from the free-quark graph of Fig. 5 is purely logarithmic. An infinite comb of in- 
finitely narrow resonances, taken with one and the same residue, see Eq. ([T0|), yields 
the same logarithm in the deep Euclidean domain. This explains why all resonance 
couphng constants are set equal in Eq. (p!0|). I hasten to add that the spectral 
density in the real world is much more contrived. The toy spectral density (|l^) is 
supposed to be a caricature, only roughly reminding the actual one. Nevertheless, 
the exercise with the toy spectral density is instructive in several respects. 



It is clear that the leading logarithmic term (|1^ carries absolutely no information 
on the mass of the lowest-lying state or on the spacing between the resonances: it 
has no built-in dimensional parameter. This information is encoded in the power 
corrections. 

Examining the expansion in Eq. ([TTl), we conclude that the high-order 

terms have factorially divergent coefficients. The factorial divergence at large n is 



due to the factorial growth of the Bernoulli numbers i?2n in Eq- (|T2|) . In QCD 
the expansion in powers of is nothing but the condensate expansion. The 

example considered teaches us that the power expansion of 11 (Q^) has zero radius 
of convergence. By itself, this is no disaster, one can work with asymptotic series 
as long as the expansion parameter is small. In our case the expansion parameter 
is r = and we are interested in the expansion in the vicinity of ~ 1. 

(Remember that we put = 1. If is substantially larger than 1, one gets very 
little information on the p meson from n(Q^).) 

In the ideal world we would calculate Il{Q^), and, hence, the spectral density 
proportional to Imn(s), exactly. In reality we must settle for n((5^), calculated in 

^In the case at hand we deal with pure powers of Our toy model is too rude to reproduce 

the logarithmic corrections and logarithmic anomalous dimensions of the condensates typical of 
QCD. 
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the Euclidean domain (positive Q"^), in the form of a truncated series. The dispersion 
relation provides a bridge between the Euclidean calculation and the spectral density, 

n(Q2) = i / + Const. (18) 

This explains why the method is referred to as the sum rules. 

The first five terms in the power expansion of the polarization operator (^) are 



InQ^ \ \ - + 

^ 2 4 3 40 



T^l/Q'. (19) 



From the structure of the series it is seen that at r ~ 1 the polarization operator 
n(Q^) is determined by its expansion with very poor accuracy, up to a factor of 2, 
at best. This accuracy is obviously insufficient to allow one to estimate, say, 
with a reasonable precision. 

One can drastically improve the accuracy by considering the Borel-transform 
of n((5^) [0]. The Borel transformation can be defined in various ways. For the 
functions obeying dispersion relations, the most convenient definition is through the 
following limiting procedure: 

i§ = lim^-A_(g2)"(^-^^ ^ g^^oo, n^oo, ^^M' fixed. (20) 

is called the Borel parameter. It is not difficult to show that applying B to 
(l/g2)" we get 

which entails, in turn 

Another trivial but useful relation we will need is 

i3(lng2) = -l. (23) 
The Borel-transformed dispersion relation (|T8|) takes the form 

n(M2) = {m{Q^)} = ^ / dslmUis) e"^/^' . (24) 

The advantages one gains in dealing with the Borel-transformed sum rules are 
obvious. First, we improve, factorially, the convergence of the power series. Thus, 
in the toy model we are playing with - an infinite comb of infinitely narrow peaks 
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presented in Eq. ([T0| ) - the l/Q"^ expansion for T{{Q^) is asymptotic, while the 
expansion of 11 (M^) is convergent] Indeed, 



,2 



N, 1 e-i/^' 
M 

and the radius of convergence of the 1/M^ expansion is 



|M;| = . (26) 

The factor 27r in the denominator ensures that we can go down to a remarkably low 
value of M^; the point ~ 1 is well inside the convergence radius. Let us examine 
the first five terms of the truncated series. 



n(M2) ^' 



127r2 



r r2 ISr^ 
1-1 ^ \ y. 

2 4 6 240 24 



r = 1 /M^ . (27) 



It is seen that in the vicinity of r ~ 1 our "theoretical prediction" is expected to 
carry an error of order ~ 2%. Such accuracy is already good enough to get the mass 
and residue of the lowest-lying state with a reasonable "resolution" . 

On the phenomenological side, proceeding to the Borel-transformed dispersion 
relations we automatically kill possible subtraction constants. What is even more 
important, the exponential weight function in Eq. (^) makes the integral over the 
imaginary part well-convergent. Thanks to the improved convergence the relative 
role of the p-meson contribution is strongly enhanced, as we will prove shortly by 
quantitative estimates. 

To see how it all works let us take a closer look at n(M^); to avoid cumbersome 
numerical factors we will change the overall normalization and will deal with 

im = ^n(M^) = ^ / dse-^/'^' Pis) . (28) 

In the toy model under consideration everything is known explicitly: the exact 
expression for /(M^), its power expansion, and the p-meson contribution to /(M^). 
The corresponding plots are displayed in Fig. 6. The exact curve has a typical shape: 
at small it is exponentially suppressed, approaching zero as exp(— m^/M^); at 
large it is flattens off and slowly approaches its asymptotic value, which is equal 
to unity thanks to a smart choice of the normalization factor in Eq. (^Sj) . This unity 
is in one-to-one correspondence with the free quark result ([T^). It is convenient to 



normalize the sum rules to the free-quark result at asymptotically large M^, and 
we will always do that. The exact curve has a steep (left) shoulder, and a shallow 
(right) one, with a maximum at slightly above m^. Remember this pattern since 
we will encounter with a very similar picture more than once in real QCD. 

The p meson mass and residue follow immediately from consideration of /(M^) 
in the small domain (the left shoulder). Indeed, in this domain all higher states 



13 



are exponentially suppressed (for instance, the contribution of the first excitation 
relative to that of p is ~ exp(— 3/M^). Therefore, /(M^) is fully saturated by p in 
the hmit 0. Fitting the left shoulder by Ci exp(— C2/M^) we would find Ci^2, 

the p-meson residue and mass, respectively. 

Alas... In the real world J(M^) is not known at small M^. All that is known 
is the large-M^ asymptotics, in the form of the expansion in the coupling constant 
q;s(M^), and a (truncated) condensate expansion. Therefore, let us pretend that in 
our toy model we have the same: a truncated power series (|27|). The corresponding 
curve is shown in Fig. 6 (where I have actually included also the O(r^) term). The 
expansion is convergent only to the right from 0.5. If we want the truncated 

series, with a few first terms included, to accurately represent the theoretical curve, 
we have to stop at 0.7 (arrow A in Fig. 6). This is the left boundary of what 

the sum-rule practitioners usually call the window or working window. 

The right boundary of the window (arrow B) is provided by the requirement that 
the first and higher excitations show up at the level not exceeding, say, 20 to 30 %. 
In this way we ensure sufficient sensitivity to the p-meson parameters. 

The larger the M^, the weaker the p-meson dominance will be. The requirement 
of the p-meson dominance forces us to move to smaller values of M^, while keeping 
control over the power expansion suggests that we move to larger values of M^. 
Thus, these two requirements are self-contradictory. A priori, in any given problem, 
it is not evident that the AB window exists at all. If it does, we can obviously use 
this fact to approximately calculate the mass and the residue of the lowest-lying 
state. The reasons why the window exists in the problems of the classical mesons 
and baryons will be discussed later. 

To get an idea of the accuracy one can expect from the sum rule calculations of 
the p-meson parameters we must investigate how sensitive our results are to various 
assumptions regarding the spectral density. The spectral density presented in Fig. 
4 is not fully realistic, for many reasons. If minor details affected the results too 
strongly, the SVZ method would be useless. The right edge of the window (arrow 
B) is chosen in such a way as to minimize the impact of the (unknown) details 
of the spectral density above the p-meson peak. Yet, it is important to obtain a 
quantitative measure of the corresponding uncertainty. 

The most obvious unrealistic feature of the spectral density (0) is the infinitely 
narrow width of all resonances from the comb. It is true that in the multicolor QCD 
the width-to-mass ratio of all hadrons built from quarks is proportional to 1/Nc 
and, thus, vanishes aX Nc ^ oo. If Nc is fixed, however, and we are interested in the 
widths of the excited states, there is another large parameter, the excitation number, 
which compensates for the effect of large N^.. Moving to the right along the s axis in 
Fig. 4, very soon we find ourselves in the situation with the overlapping resonances 
- the 5 functions in the comb representing the spectral density are smeared and the 
spectral density becomes perfectly smooth. Let us discuss the consequences of such 
dynamical smearing. 

First of all, we must establish the dependence of the width r„ on the excitation 
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Figure 6: The sum rule for the p meson in the toy model. 



number n. If one believes that a string-like picture of color confinement develops in 
QCD, an asymptotic estimate of r„ at large n becomes a simple exercise. Unlike the 
p meson and other low-lying classical states, for highly excited states the string-like 
picture of the chromoelectric tubes is expected to be fully relevant. 

When a highly excited meson state is created by a local source, it can be consid- 
ered, quasiclassically, as a pair of (almost free) ultrarelativistic quarks; each of them 
with energy m„/2. These quarks are created at the origin, and then fly back-to- 
back, creating behind them a flux tube of the chromoelectric field. The length of the 
tube L ~ rrinl where represents the string tension. The decay probability is 
determined, to order by the probability of producing an extra quark- ant iquark 
pair. Since the pair creation can happen anywhere inside the flux tube, it is natural 
to expect that 

r„, ~ i-LA^ = (29) 

where i? is a dimensionless coefficient of order one. I used here the fact that within 
the quasi classical picture the meson mass m„ ~ LIS?' . Thus, the width of the n-th 
excited state is proportional to its mass which, in turn, is proportional to ^Jn for 
the linear Regge trajectories Q 

This simple estimate was obtained in Ref. long ago. Since the argument 



bears a very general nature, the square root dependence of r„ should take place in 



^ Let us note in passing that the l/iV^ corrections due to creation of two quark pairs are of 
order L? jN"^ within this picture. Since L ~ m„ ~ ^Jn, the expansion parameter is ^/njN^. 
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all models with linear confinement. And it does, indeed! Recently, the square root 
formula for r„ was obtained numerically in the 't Hooft model [jl5 . 



If Eq. ( P^ ) does indeed take place, it is not difficult to find the impact of the 
resonance widths on the spectral density. The infinitely narrow pole is substituted 
with 

1 1 1 



1 1 



(30) 



where 



and I used the fact that 1/Nc is a small parameter. This explains the last transition 
in Eq. (^). For the same reason, in the third transition I replaced m\/Nc with 
N~^Tt'^Q'^\nQ'^ . Near the pole (i.e. at = —mfj both expressions give the same 
in the leading l/N^. approximation, and are equally legitimate 0. 

As a result, all poles, which lied previously at positive real g^, are now shifted 
onto unphysical sheets, away from the physical cut (Fig. 7). The 6 functions in the 
spectral density are replaced by peaks with finite widths. Correspondingly, the toy 
model number two, which replaces Eq. (0) is 

n(Q') = H^) + Const. , (32) 

where now z is given by the following formula 



- . (33) 

The additional factor (1 — 7)^^ in the overall normalization is chosen to reproduce 
the free quark result (|np at large Euclidean Q"^. Note that the polarization operator 
defined by Eqs. (|3^ and (^) has the correct analytical structure: it is non-singular 
in the whole complex plane, apart from the cut at real positive g^. The singu- 
larities associated with the poles of the ip function lie on the unphysical sheet (Fig. 
7). 

The spectral density (i.e. the imaginary part of 11 at = s + ie) takes the form 
of the sum of (modified) Breit-Wigner peaks 

( \ _ 3 ^ s^~^sin7r7 
^^'^^ ~ 7r(l -7) ^0 [s^-T cos7r7 - (1 + 3A;)]2 + [s^-t sin n-f]'^ ' ^ ^ 

It is depicted in Fig. 8 for i? = 0.6 and = 3. Although the choice of the 
parameter B is rather arbitrary, it is not unreasonable; it corresponds to F/m ~ 0.2. 
Experimentally the p meson has a close width-to-mass ratio. 

^The logarithm introduces not only the imaginary part, i.e. the width, but a shift in the real 
part, equivalent to a shift in the resonance mass. Both effects are of order l/N^- 
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Upper side of the cut 



(a) (b) 

Figure 7: Analytical structure of the polarization operator, (a) The polarization 
operator must be analytic everywhere in the complex plane, except the cut 
running on the negative real semi-axis of (positive real semi-axis of q^). The 
imaginary part of the polarization operator must be positive at the upper side of 
the cut; {b) The mapping of the plane onto the z plane, Eq. (^). The physical 
sheet on the plane corresponds to the z plane with the shaded sector removed. 
The boundaries of the sector correspond to the lower and upper sides of the cut. 
The angle of the removed sector is proportional to 7. 



17 



P(s) 



toy model # 2 
(modified) Breit-Wigner 




Figure 8: The spectral density corresponding to a sum of (modified) Break- Wigner 
peaks with one and the same width-to-mass ratio 0.2. All residues are set equal, 
and 7 ~ 0.06369, see Eq. (|3lD. The s axis is in the units m^. 



We can see on the picture a pronounced p-meson peak, accompanied by a deep 
dip, and then an almost smooth curve which approaches the asymptotic value (unity) 
in an oscillating mode. The smearing of the 6 functions into a smooth curve occurs 
because the width of the excited states is proportional to m„ while the distance 
between two neighboring states Am„ ~ l/'^^n- It is assumed, of course, that A*",; 
is fixed. The larger the value of Nc, the further we must go, to higher excitation 
numbers, for the resonances to overlap. Figure 8 shows that at A^^^ = 3 a smooth 
curve starts right after the ground state. The first excitation already belongs to 
the smooth curve. Remember this pattern of the spectral density - a conspicuous 
first peak accompanied by a smooth (and oscillating) curve approaching its constant 
asymptotic value. The gross features of this picture follow from very general theo- 



retical arguments ||T6|. We will see shortly that experimental data, being different 
in fine details, exhibits the very same type of behavior. 

The spectral densities shown in Figs. 4 and 8 at first sight do not produce 
an impression of close relatives. Let us integrate them over, with the exponential 
weight function, and compare J(M^) obtained in this way. The comparison is shown 
in Fig. 9. The dotted curve corresponds to the comb of S functions. It was already 
presented in Fig. 6. The solid curve corresponds to Eq. (|3^ and Fig. 8. Finally, 
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Figure 9: J(M^) for three models of the spectral density. 

the dashed curve is obtained with an extremely crude model of the spectral density 
depicted in Fig. 10, 

p(s) = 35(s - 1) + e{s - So) , So = 2.7 . (35) 

It presents an infinitely narrow p, plus a gap, plus all higher excitations fused into 
"continuum" which starts abruptly at s = Sq and coincides (at s > Sq) with the 
asymptotic free-quark expression for the spectral density. This is the original SVZ 
model [||]. It, obviously, captures only gross features, and misses all details: the fact 
that the gap at s > is not quite empty, the onset of continuum is not so abrupt, 
and the limit p{s) — 1 at s — » oo is achieved through oscillations. Nevertheless, 
it is seen that I{M'^) in all three models come out very close to each other. The 
Breit-Wigner model of Eq. (p4D gives a little bit less steep fall off to the left of the 
window, at —>■ 0. This is quite understandable, since in this model there is a 
non- vanishing spectral density at s < m^, due to the p meson width. 

Comparison of three curves in Fig. 9 teaches us a lesson: by examining I{M'^), 
calculated approximately inside the window, one cannot expect to predict, with any 
reasonable accuracy, such fine structure of the spectral density as the width of the 
p meson, and the shape of the curve in the continuum (i.e. higher excitations). 
As Fig. 9 shows, it is hardly possible to distinguish between three models under 
discussion, which, in a sense, represent extreme situations. The best one can hope 
for is estimating m^, its residue, and, to a lesser extent, an effective value of s where 
the dip following the p-meson peak ends (i.e. sq). 
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Figure 10: The original SVZ model of the spectral density. 



In a broader context, this message is instructive for the lattice practitioners too. 
Doing calculations in the Euclidean domain - a common feature of the sum rules 
and lattice technology - derives virtually no information about the spectral density 
beyond the ground state in the given channel. Indeed, the original SVZ model and 
the comb of the 6 functions, representing drastically different p{s) point-by-point, 
lead to I{M'^) coinciding with each other to a high accuracy. To distinguish between 
these two scenarios one needs exponential precision going far beyond the level so far 
achieved in the lattice calculations. 

Now we are ready to explore the sensitivity of /(M^) to the p-meson parameters. 
To this end we will take the SVZ model in the form 



p(s) = Ci6{s - C2) + 9{s - So) (36) 

and let the parameters float by, say, 10% around their "reference" values, Ci = 
3, C2 = 1 and sq = 2.7. We then compare I{M'^) obtained in this way with I{M'^) 
emerging in our (most realistic) Breit-Wigner toy model, see Eq. (|3^) . 

Figure 11 illustrates the sensitivity of the sum rule to the p meson residue (cou- 
pling constant). The upper and lower dashed curves correspond to Ci = 3.3 and 2.7, 
respectively. Figure 12 illustrates the sensitivity to the p meson mass. The upper 
and lower dashed curves correspond to = 0.9 and 1.1 of the experimental value, 
respectively. It is seen that it is quite realistic to expect to get g"^ and with 
the accuracy ~ 10% by inspecting li^M"^) in the window and fitting the theoretical 
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Figure 11: Changing the p meson residue. The sohd curve corresponds to a "refer- 
ence" spectral density presented in Eq. (plf). The dashed curves correspond to Eq. 
(^) with C\ = 3.3 and 2.7, respectively. The reference value is Ci = 3. 

prediction by the model (|3^) . The estimates of sq are less accurate. One should not 
be surprised, of course, since the sum rules were designed to be most sensitive to 
the ground-state parameters and relatively insensitive to the details of the spectral 
density in the continuum. The exponential weight in the definition of /(M^) takes 
care of this feature. Deviations of sq from its reference value at the level of 10% 
lead to rather insignificant changes in J(M2). Deviations of sq at the level of 20% 
are quite noticeable. 

I hasten to add, however, that the 10% accuracy is not a typical outcome of the 
sum rule analysis. The p-meson channel is most favorable from the point of view of 
applications of the SVZ sum rules. In a sense, this is a dream case: the role of con- 
tinuum with respect to p is as tempered as it can possibly be, and higher (unknown) 
condensates in the truncated condensate expansion show up at remarkably low val- 
ues of so that the working window is comfortably wide. In many other channels, 
mesonic and baryonic, we have to deal with a narrower window. The general rule 
is: the narrower the window the worse the accuracy. In some channels, as we will 
see later, the window shrinks to zero. Then the SVZ method fails. The reasons why 
it is successful in some cases and fails in others, when properly understood, give us 
a unique hint as to the structure of the QCD vacuum (see Sect. 8). 
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Figure 12: Changing the p meson mass. The dashed curves correspond to the 10% 
variation of m^, around its reference value (m^ = 1). 



3 Vacuum Condensates 

It is high time to abandon our baby version of QCD and proceed to the real thing. 
What can be said about /(M^) now? 

Needless to say that no analytic calculation of J(M^) at < A^, that would 
be based entirely on the first principles, exists. QCD at large distances is not 
yet solved analytically. Under the circumstances it seems reasonable to start at 
short distances and advance to larger distances "step by step" , staying on the solid 
ground - working with the microscopic variables, quarks and gluons, where it is fully 
legitimate. At short distances the quark-gluon interactions are adequately described 
by perturbation theory. Certainly, as we have just learned by inspecting the toy 
models, the (truncated) perturbative series is not going to yield us any estimates 
relevant to the p meson parameters. To get them we should add at least some 
information regarding the large distance dynamics. 

In the ideal world this information would be obtained from the theory per se. 
In the real world we may try to parametrize the effects caused by the vacuum 
fields. If the quark-antiquark pair injected in the vacuum by the current does 
not propagate too far, its impact on the vacuum fields is, hopefully, not drastic. This 
means that the polarization operator can be well approximated by the interaction 
of the valence quarks with a few vacuum condensates. For the validity of this 
assumption it is necessary that the characteristic frequencies of the valence quarks 
inside the p meson u be larger than the characteristic scale parameter of the vacuum 
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medium /x. Certainly, we can count only on an interplay of numbers: since the only 
dimensional parameter of QCD is A, all quantities of dimension of mass are of order 
A. It may well happen, however, that, say ~ 3A while fi ~ A/2. In a sense, the 
success of the SVZ sum rules confirms this assumption a posteriori. 

On general grounds it is expected that all gauge invariant Lorentz singlet local 
operators built of the quark and/or gluon fields develop non- vanishing vacuum ex- 
pectation values (VEV's). It is convenient to order all operators (and their VEV's) 
according to their normal dimensions. The higher the dimension, the higher the 
power of 1/M of the corresponding coefficient. This fact allows one to control the 
condensate expansion inside the working window. 

The lowest dimension (zero) belongs to the unit operator. The unit operator is 
trivial. When one calculates the perturbative contribution to /(M^), one actually 
calculates the coefficient of the unit operator. 

The operator of the next lowest dimension (three) is the quark density operator, 

0, = qq. (37) 

No other gauge and Lorentz invariant operators of dimension three exist. 
There is one operator of dimension four, the gluon operator 

= —G'^^u^% ■ (38) 



The factor a^/Tr appears in the definition naturally. Since both operators, (|37|) and 
(p8|), play a very special role we will interrupt our excursion before passing to higher 
dimensions, to make several important remarks. 

In the chiral limit (i.e. when the quark mass term in the Lagrangian is put to 
zero) Oq is the order parameter - its VEV signals the spontaneous breaking of the 
axial SU{Nf) symmetry and the occurrence of the corresponding massless pions. In 
perturbation theory the chiral symmetry remains unbroken, of course, (Og) vanishes 
identically. It is tempting to say then that (Og) ^ measures deviations from 
perturbation theory, and the same refers to all other condensates. 

In the early days of non-perturbative QCD such an understanding was widely 
spread. I hasten to say that this is a wrong understanding. It is impossible to define 
the condensates as "truly non-perturbative" residues obtained after subtracting from 
them a "perturbative part" [[l^]. A heated debate took place in the literature in the 
eighties regarding the possibility of defining the condensates, in a rigorous way, by 
subtracting a "perturbative part" . Glimpses of the debate can be seen in Ref. flEj 
where it was clearly emphasized, for the first time, that the procedure of isolating 
"purely non-perturbative" quantities does not (and must not) work 0. I will outline 

^Surprisingly, the debate was revived recently, over a decade after the issue had been seemingly 
settled, in connection with the heavy quark theory. This theory, as it exists today, is also based 
on the operator product expansion and is a close relative p9[ | of the SVZ sum rules. One can ask 
there the very same questions that are usually raised in connection with the condensate expansion 
in the SVZ method. Later on, in Sect. 10.2, we will return to the issue and discuss some elements 
of the heavy quark theory. 
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the proper procedure below using the gluon condensate as the most instructive 
example. 

Unlike Og, the gluon operator Oq is not an order parameter. There is no known 
symmetry whose spontaneous breaking would generate a vacuum expectation value 
(Og)- Correspondingly, (Og) is generated in perturbation theory. Physically, (Og) 
measures the vacuum energy density Svac- To see that this is indeed the case let us 
use the fact that the trace of the energy-momentum tensor 

(^'il = ■ (39) 



This expression is nothing but the scale anomaly of QCD [^. Strictly speaking, the 



coefficient on the right-hand side contains higher orders in a^; they are not essential 
for our purposes, however, and will be ignored. The vacuum expectation value of 
the left-hand side is obviously 4£vac- Hence 

evac = (Og) . (40) 

Everybody knows that the vacuum energy density badly diverges in field theory, 
as the fourth power of the cut-off 0. If so, how can one make sense of the gluon 
condensate? 

The vacuum fields fiuctuate; these fluctuations contribute to the vacuum energy 
density. High frequency modes of the fluctuating flelds belong to the weak coupling 
regime; their contribution in the correlation functions is described by perturbation 
theory, giving rise to the standard perturbative expansions. What we are interested 
in are the low frequency modes, soft vacuum flelds that are responsible for the 
pecuhar properties of the vacuum medium. By saturating the condensate by the 
soft modes only we get a consistent deflnition of the condensates that automatically 
solves the problem of divergencies. 

Speciflcally, we must introduce a somewhat artiflcial boundary, all fluctuations 
with frequencies higher than /i are supposed to be hard, those with frequencies lower 
than /i soft. The parameter /i is usually referred to as normalization point. Only 
soft modes are to be retained in the condensates. Thus, the separation principle 
is "soft versus hard" rather than "perturbative versus non-perturbative" . This is 
the basic principle of Wilson's operator product expansion, which, in turn, is the 
foundation of the SVZ sum rules. Being deflned in this way the condensates are 
explicitly /i dependent. All physical quantities are certainly fi independent; the 
normalization point dependence of the condensates is compensated by that of the 
coefficient functions. 

Generalities of Wilson's approach and its particular implementation in QCD will 
be further discussed in Sects. 5 and 6. Here we complete our consideration of the 
gluon condensate. 



^ Let me parenthetically note that in supersymmetric glyodynamics the vacuum energy density 
vanishes, and so does the gluon condensate. 
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Figure 13: The one-gluon correction to the D function. 



It is instructive to trace how the gluon condensate emerges from a shghtly dif- 
ferent perspective. A natural starting point is the diagram of Fig. 13 presenting the 
one-gluon correction to Il{Q^) (two other similar graphs with the one-gluon exchange 
are not shown). For convenience we will deal with the so called D function, 

D{Q^) = -{47r^)Q\dU/dQ^) , (41) 

rather than with the polarization operator 11 itself. In terms of D{Q^) the bare 
quark loop of Fig. 5 becomes just unity. 

If the gluon momentum is denoted by k the one-gluon term in D[Q^) takes the 
form 

AD{Q^) =a,j edk' ^ F{k\ (42) 

where it is implied that the angular integration over the orientations of the vector 
k with respect to q is already carried out. The factor 1/fc^ comes from the gluon 
propagator, while Fik"^, Q^) represents the rest of the graph. The function F{k'^, Q"^) 
was calculated in Ref . |^ (it can also be extracted from Appendix B of Ref . ) , 



F{k\ Q') = {(i - 1^^) ^ + (1 + ^) iL2{-r) + lnrln(l + r)]} , (43) 



3n Q 

where 
and 



rx 

L2ix) = - dyy-'^\n{l - y) 
Jo 



is the dilogarithm function. The analytic continuation of the function (^31) to r > 1 
can be more conveniently rewritten as 

(1 + r) [L2(-l/r) -lnrln(l + r-i)]} . (44) 

Substituting Eqs. (|^) and (|4^ ) in Eq. ( ^2]) and doing the integral over k"^ in 
the most straightforward manner one gets AD = ag/ir, the standard well-known 
one-gluon correction to the D function. 
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Figure 14: The distribution function Q F{k\Q ) defined in Eq. p^) versus r. 



The plot of the function F{k'^^ Q"^) is shown in Fig. 14. Although it has a clear- 
cut peak at r ~ 0.5, it presents a distribution spanning the entire range from = 



to k = oo. The gauge coupling constant as in Eq. 



is literally a constant 



provided one limits oneself to the graph of Fig. 13. In higher orders, however, ag 
starts running, and it is evident that — > Then as{k'^) should be placed 

inside the integral in Eq. (^) rather than outside. By doing so we effectively resum 
a subset of the multi-gluon graphs. 
Certainly, putting 



1 + 



Air 



1 



(45) 



inside the integral in Eq. (^2]) we do not account for all higher-order corrections. 
Only those are included that are responsible for the running of as in the leading 
logarithmic approximation. Usually they say that the corresponding multigluon 
graphs are of the bubble chain type (Fig. 15). Strictly speaking, in the covariant 
gauges the gauge coupling renormalization ( |i5| ) in QCD (unlike QED) does not 
reduce to the bubble insertions in the gluon propagator depicted in Fig. 15. Various 
graphs of a more complicated structure are involved too; they combine together to 
produce a gauge- invariant expression (^5]). The bubble chain saturation of Eq. (^5]) 
takes place in the physical gauges, e.g. in the Coulomb gauge. In the covariant 
gauges one can artificially reduce the problem to the bubble chains, without explicit 
identification of the full set of graphs, by exploiting the so called "large negative 
Nf trick". We will not dwell on this technical issue, since it is irrelevant for our 
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purposes. The interested reader is referred to Refs. 

Even though the bubble chain is a very specific subset of graphs, and the majority 
of the multigluon graphs are left aside in this procedure, it still makes sense to 
substitute the running as(fc^) inside the integral in Eq. (|42|), since in this way we take 



26 



into account an essential part of the underlying dynamics: the fact that the quark- 
gluon interaction becomes weak if the gluon is far off-shell, and, on the contrary, 
becomes stronger as the gluon virtuality decreases. If we follow this strategy, then 



dk^ 



1 + 



47r 



In- 



F{k\Q^). 



(46) 



We immediately run into a problem here. The running as has the Landau pole at 
small /c^, the square bracket in the integrand explodes right inside the integration 
interval. The explosion occurs at 



k"^ = exp 



A' 



(47) 



Although the function F{k'^,Q'^) is suppressed at small fc^, i.e. r — > 0, (see Fig. 
14), it by no means vanishes at r = A^/Q^, rendering the integration impossible. 
In order to do the k'^ integral literally, as it is given in Eq. (^), we must say how 
the singularity at k'^ = A^ is to be treated. Some people say: "let us take the 
principal value", others bypass the singularity by shifting the integration contour 
in the complex k"^ plane, still others do something else. All these prescriptions are 
arbitrary and physically meaningless for obvious reasons. Indeed, let us introduce 
the normalization point fi = several units xA and split the integration over k"^ in 
two parts: from zero to /i^ in the first integral and from /x^ to oo in the second. 
Then at k'^ > fi"^, in the weak coupling domain, we do trust Eq. (^) while below /i^ 
we do not. At best, below //^ one can use Eq. (^61) for the purpose of orientation. 
At small fc^ 



(48) 



implying that 



TT 



^ k'dk' 



1 + 



Air 



1 



1 

b 



^ ^ k'de 



A^ 



Stt A^ 



P-A2 ^ ^^^^ 

This is a very crude estimate; I substituted the square bracket by the corre- 
sponding pole, and then evaluated the integral by equating it to its imaginary part, 
i.e. substituting (P - A^)"! nS{k'^ - A^). 

From this exercise we learn the following lesson. The domain fc^ < fi"^ gives a 
contribution to D{Q'^) which scales as l/Q^. The coefficient in front of cannot 
be reliably calculated and must be parametrized by the gluon condensate - that is 
the best we can do. It is not accidental that the /c^ expansion of F{k'^,Q'^) starts 
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Figure 15: One of the graphs from the bubble chain set. 



from k'^, while the term of the zero-th order is absent. If it were not the case, the 
condensate series would have to start from 0{Q~^). This is impossible, however, 
since no gauge invariant local operator of dimension two exists in QCD. The absence 
of such operator and the absence of terms in the distribution function F{k'^, Q^) 
are in one-to-one correspondence. 

At first sight, one could try to circumvent the problem of the Landau pole inside 
the integration domain by expanding the square bracket in Eq. (^) in the series 
in as{Q^)- This is one of the "remedies" sometimes cited in this context. I put the 
word remedy in the quotation marks since in fact it does not work. True, in every 
given order of the expansion the integral P6| ) becomes well-defined. It is easy to 
see, however, that the expansion is going to be factorially divergent in high orders. 
Indeed, 

fcWfln^) ~2-"n!, n > In ^ . (50) 
\ k^ J 

It is seen that perturbation theory per se carries the seeds of the gluon con- 
densate. The perturbative series cannot be consistently defined unless the gluon 
condensate is introduced. Once we cut off the k^ < /i^ tail of the integral from 
the perturbative expression obtained from Eq. (^Bf) , the factorial divergence of the 
as series at high n disappears, the series becomes well-defined, and so is the gluon 
condensate. 

The factorial divergence of the perturbative series of the type we have just dis- 
cussed is called the infrared renormalon ^ for the reasons which need not 
concern us here. It is associated with the bubble chain graphs of Fig. 15 and is 
inevitable if one forces the perturbative integrals to run all the way down to A;^ = 0. 
Bounding the integration interval in the perturbative part from below by yU^ we get 
rid of the infrared renormalon altogether. The domain below /i^ is not lost; it is 
fully represented by the gluon condensates. Higher condensates appear as higher 
order terms in the fc^ expansion of F{k'^,Q'^). The intimate relationship between 
the infrared renormalons and condensates was first revealed by Mueller p5| . 



I hasten to warn that one should not literally equate the infrared renormalons 
to the condensate expansion. Even if one sticks to a certain particular definition 
for summation of the bubble chains, that eliminates ambiguities, say the principal 
value prescription in the integral (|46|), one typically gets a numerical value of the 
term which is grossly off compared to the gluon condensate, let alone ambigu- 
ities associated with various possible choices of the regularization prescription. For 
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instance, in the case at hand, 

^(Q') = 1 + ^((^g) + ..., (51) 

and the 0{1/Q'^) term is 0.08 GeV^/Q'^ under the standard choice of the gluon con- 
densate (see below). Equation (|4^) yields, on the other hand, only 0.004 GeV^/Q^ at 
A = 0.2 GeV (0.01 GeVVQ^ at A = 0.25 GeV). Within the SVZ method, based on 
Wilson's expansion, the role of the infrared renormalons is purely illustrative. They 
show, in a very straightforward manner, that without the condensates the pertur- 
bative series cannot be defined and, on the contrary, introducing the condensates 
one eliminates ambiguities associated with the high-order tails of the perturbative 
series []. 

On the other hand, in some (perturbatively infrared stable) processes that have 
no operator product expansion, most notably in the jet physics, the infrared renor- 
malons, in spite of all their obvious limitations, become poor man's substitute for 
the condensate expansion. Although, unlike OPE, the renormalon analysis does 
not provide us with the proper numerical coefficients, we may infer from the renor- 
malons what powers of 1/E (or 1/Q) appear in the nonperturbative parts of such 
quantities as, say, thrust. Whether the nonperturbative corrections are 0{1/E) or 
0{1/E'^) is clearly a question of paramount importance for phenomenology. Using 
the infrared renormalons for counting the powers of 1/E in the processes without 
OPE is a totally fresh idea which can be considered as a distant spin off of the SVZ 
method, see Ref. for a review. Previously the nonperturbative corrections were 



just ignored since nobody knew how to approach the issue scientifically. 

Now, when the meaning of the condensates appearing in Wilson's expansion is 
hopefully clear, we can proceed to a brief discussion of their numerical values. In 
principle, two alternative approaches are conceivable: (i) calculating the condensates 
from first principles (e.g. on the lattices); (ii) extracting the condensates from the 
sum rules themselves. Attempts of the lattice calculation of the condensates were 
reported in the literature (for a review see [^). The main difficulty, which, to a 
large extent, devaluates this work is the fact that the strategy of "isolating a non- 
perturbative residue from the perturbative background" was adopted. As we already 
know, this strategy is doomed to failure. Instead, one should have consistently 
implemented Wilson's procedure. This has never been done, however. Some initial 
ideas as to how Wilson's procedure can be implemented on the lattices are presented 



in Ref. [||]. 



^The infrared renormalon is not the only source of the factorial divergence of the as series in 
QCD. The expansion of Eq. ( ^ ) contains also the so called ultraviolet renormalon - factorial 
divergence of high-order ag terms coming from the large fc^ domain. This factorial divergence is 
readily summable, however, as is perfectly clear from the unexpanded expression (p6[). Indeed, the 
integrand has no singularities at large and is well convergent. Below, the ultraviolet renormalons 
will never be mentioned again. Those readers who want to familiarize themselves with this subject 
are referred to an excellent review p6| . A factorial divergence of a different nature showing up in 
the coefficient functions is briefly discussed in Sect. 5. 
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Therefore at present, as twenty years ago, one has to exploit the sum rule them- 
selves in order to determine the condensates. One picks up certain channels where 
sufficient experimental information on the corresponding spectral densities is avail- 
able. These channels are sacrificed, i.e. the analysis goes in the direction opposite 
to conventional: the values of the condensates are extracted from the Euclidean 
correlation functions obtained through dispersion integrals. The ffist determination 
of the gluon condensate was carried out exactly in this way, from the sum rules in 
the J/tp channel [^, 0]. The corresponding value will be referred to as standard, 

{Og) = 0.012 GeV^ . 

Since then this estimate was subject to multiple tests, sometimes with conflicting 
results. A brief discussion can be found in the beginning of Sect. 2 in Ref. 0. 
The overall conclusion is that certain deviations from the standard value (say, at 
the level of ~ 30%) are not ruled out 0. One of the reasons limiting the accuracy 
is the fact that the gluon condensate was determined within a simplified (the so 
called practical) version of the operator product expansion, which is strictly speaking 
somewhat ambiguous. We will dwell on the issue in Sect. 6. At the present level 
of understanding it is highly desirable to repeat the analysis within the framework 
of consistent Wilson's procedure. This is not done so far. To take into account a 
certain degree of numerical uncertainty it is reasonable to accept that 

{Og) = 0.012 xg GeV\ (52) 

where a dimensionless numerical factor xg is allowed to float in the vicinity of unity. 
This parametrization will be used below in the p-meson sum rule. 

Let us pass now to the quark condensate (Og). As was mentioned, this con- 
densate is the order parameter for the spontaneous breaking of the chiral symme- 
try. Hence, its value can be independently determined from the corresponding phe- 
nomenology. In particular, the celebrated Gell-Mann-Oakes-Renner formula relates 
the product of the light quark masses and {Og) to observable quantities, 

(m„ + ma) {uu + dd) = -Mlfl , (53) 

where is the pion constant (/^r ~ 133 MeV) and M^^ is the pion mass. The original 
estimate of {Og) used in Ref. was obtained by substituting the crude estimates 
of the light quark masses that existed at that time [^, |3! 



In this way one arrives at 



mu + irid^ 11 MeV . 

(Cg) = -(250 MeV)^ (54) 



®I leave aside extremist and to my mind unfounded statements in the literature, that the 
standard value underestimates the gluon condensate by a factor of 2 to 5. 
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Since then considerable work has been carried out to narrow down the theoretical 
uncertainty. First of all, chiral corrections to the Gell-Mann-Oakes-Renner formula 
were derived and analyzed. Moreover, much effort has been invested in perfecting 
our knowledge of the quark masses. The progress is summarized in Ref . , where 
an extensive list of references can be found. No dramatic changes occurred. The 
value of rUu + ma extracted from phenomenology of the chiral symmetry breaking 
went up by about 30%, with an error at the level of 15%. Thus, one is tempted to 
say that the actual value of the quark condensate is close to its "standard" value 
quoted in Eq. (^); perhaps, somewhat suppressed |. 



Events took a dramatic turn in 1996, after results of the lattice calculations of 
the light quark masses (with two dynamical quarks) became known. The lattice 
result for m„ + lies significantly lower (by a factor of two to three, see e.g. ^^^) 
than any of the reasonable analytic estimates! Were these results confirmed, this 
would mean that the quark condensate is significantly larger than the number given 
in Eq. (H). 

I would caution against hasty conclusions at this point. The discrepancy between 
the analytic and lattice estimates may well be an artifact of the lattice procedure. 
Putting dynamical almost massless quarks on the lattice is a notoriously difficult 
task. On the other hand, the discrepancy, if persists, may prove to be a signature of 
a very interesting physical phenomenon - strong dependence of physical quantities 
on the number of light dynamical quarks. The surprisingly low lattice result for 
mu + md was obtained under the assumption that the number of light quarks is two, 
while in fact it is three. We will return to the issue of possible strong dependence 
on Nf in Sect. 12. 

Under the circumstances it seems reasonable, analogously to the gluon conden- 
sate, to parametrize the quark condensate as 

(C,) = -(250 MeV)^Xg, (55) 

and let Xq vary around unity. 

The full catalog of all relevant operators up to dimension six, was worked out 
in Ref. 0]. It is not as large as one might think at first sight. There is only one 
operator of dimension five, the mixed quark-gluon operator 

OqG = igqf^fivG^yq (56) 

where G^i, = G'^^^T'^ (T" are the color generators) . This operator plays no role in the 
p-meson sum rules, to be considered below: because of its "wrong" chirality it can 

^The quark mass is definition dependent. Usually the MS scheme and a "reference" normal- 
ization point n — 1 GeV are implied. Not to confuse the reader we will follow this convention, in 
spite of shortcomings of the MS scheme. One cannot avoid explicitly introducing the normaliza- 
tion point /i in the case of the quark condensate: Oq has a non-vanishing anomalous (logarithmic) 
dimension. The qq vertex dressed by gluons is logarithmically suppressed. Including the gluons 
with off-shellness from Mq down to /i suppresses the vertex by a factor cx [as{n)/as{Mo)]'^/^ . This 
formula explicitly demonstrates why /i cannot be put to zero in the case at hand: the logarithm 
explodes. 
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enter only being multiplied by the quark mass m^. It is very important, however, 
in a wide range of problems involving baryons |^ , and mesons built from one light 
and one heavy quarks [^. 

At the level of dimension six, there is one operator built from three gluon field 
strength tensors (symbolically GGG) and several four-quark operators of the type 

Oi, = iqiT,q2){q3T2q^) (57) 

where Fi 2 denote certain combinations of the Lorentz and color matrices and 
{i = 1,...,4) are the light quark fields of different flavors, u,d and s. The overall 
flavor must be conserved of course, but qi need not coincide with q2 and so on. 

The three-gluon operator is expected to have a significant impact in heavy 
quarkonium ; it does not appear in the p- meson sum rules, as was first shown in 



Ref. I^H], by virtue of a very elegant theorem. 

As for the four-quark operators, their vacuum matrix elements are not known 
independently. It became standard to evaluate them in the factorization approx- 
imation, which was first applied in this context in Note that using the Fierz 
identities for the color and Lorentz matrices one can always arrange the operators 
(^4^ to have a natural color flow. By natural I mean that Fi 2 do not contain color 
matrices, and the color indices of quarks are contracted inside each bracket in Eq. 
(^). The color of qi flows to q2 and that of q^ flows to q^. If the number of colors 
Nc were large, then the expectation value of O^g would be totally saturated by the 
vacuum intermediate state, 

= (giFig2)(g3r2g4). (58) 

Corrections to this equation are formally of order of It is clear that in the 

factorization approximation, after the natural color flow is achieved by the Fierz 
rearrangement, the only surviving structure is that with Fi^2 = 1- 

In the real world Nc = 3, and we certainly do expect deviations from factorization 
at a certain level. As far as we can tell today, factorization works surprisingly 
well at least for those four-quark operators that appear in the p-meson sum rule. 
All attempts to detect deviations, both on the lattices ||39l and in the sum rules 



themselves |^0|], gave results consistent with zero, within errors that typically lie 
in the 10% ballpark. This is quite nontrivial, considering the fact that quite a 
few examples are known where large numerical coefficients neutralize formal 
suppression factors. Whatever the reasons might be, we will accept Eq. (|58D in 
what follows. 



4 p Meson in QCD 

Experience accumulated in Sect. 2 will now be applied for explorations in actual 
QCD. As we learned from the toy models, the polarization operator n(Q^) (or its 
Borel transform n(M^)) in the deep Euclidean domain has an expansion in InQ^ 
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(perturbation theory) plus power terms or (l/M^)*^ (the condensates). If in 

the toy models the power terms (truncated in certain order) can be made as large or 
as small as we want, in QCD they are determined dynamically; through interactions 
of the valence quark pair with the soft vacuum fields. This interaction gives rise 
to the condensate expansion. The first power term is due to the gluon condensate, 
while the next one is associated with the four-quark condensate. In principle, higher 
condensates were analyzed in the literature too, but we will not discuss them in this 
Lecture. 

The quark-antiquark pair qq with the total energy ^/s is injected in the QCD 
vacuum either by the virtual photon {e^e~ annihilation) or W boson (r decays). 
The qq pair, being injected, starts evolving according to the dynamical laws of QCD. 
At first the quarks do not feel the impact of the vacuum "medium" . As separation 
between them grows, the effects of the medium become more and more important, 
so that eventually it prevents quarks from appearing in the detectors. The injected 
quarks get dressed and materialize themselves in the form of hadrons. In the sum 
rule approach we control only the beginning of this process. 

In the sum rule framework, the condensates become important in the vicinity of 
the left edge of the window. To the left of the arrow A they explode, and theoretical 
control is lost. On the other hand, the value of I{M'^) near the right edge of the win- 
dow is determined mostly by ordinary perturbation theory. All perturbation theory 
sits in the unit operator. Technically it is convenient to formulate the corresponding 
calculation in terms of the perturbative corrections to the spectral density. We take 
the graph of Fig. 5, and attach to it various gluon and quark loops (e.g. Figs. 13 
and 15). In this way we get the perturbative part of p(s); at present in the vector 
isovector channel it is known up to third order in as{s) [^. In the MS scheme the 
result takes the form 

p(s) = 1 + A;^a + k2a^ + ha!^ (59) 

where a{s) is given in Eq. (|^). In the low-energy domain the number of active 
flavors is Nf = 3. For Nf = Nc = 3 the parameter 6 = 9, and the values of ki are 

4 

ki = k2 = 0.729A;i , fcg = -2.03fci . (60) 
9 

Equation (|59D can be immediately translated into the prediction for the per- 
turbative part of /(M^). The perturbative expansions for J(M^) and p{s) do not 
coincide. To get the former one must integrate the latter with the exponential weight 
function. The transition is readily carried out with the aid of the formula 



ds exp(— s/M^ 



/ 



1 



A/2 



(61) 

where 7 denotes Euler's constant, 7 = 0.577. Expanding this expression near integer 
values of u one also gets necessary expressions for the integrals containing logarithm 
of logarithm (In In). 
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In this way we arrive at 



/pc(M^) = l + kia 



eT 



M2 

eT 



n 2 



TT 

"6 



M2 

e'' 



n 3 



(62) 



Similar expressions were obtained in Ref. were a(M^/e'^) was reexpressed in 

terms oi a{M^). It is more convenient to work directly with the expansion (|62|). 

The quark mass corrections in p{s) and I{M'^) are also known, they are propor- 
tional to rf/M^ and are negligibly small due to the smallness of the current quark 
masses. They will be discarded since the corresponding uncertainty is invisible in 
the background of other uncertainties. 

Now we return to the condensates. The modern perturbative calculations of the 
coefficients C„ are based on the background field technique. This is an important 
aspect of the SVZ approach as it exists today. The method is quite elegant, to say 
nothing that it is much more economic than the original "brute force" calculations of 
Ref. [m . We will not discuss the technique per se, however. Conceptually everything 
is clear here. If you need to do a calculation, you just go and learn the corresponding 



technique using, for instance, the review paper |44 



The coefficients Cg and C^q, to the leading order, are known from the ancient 
times fBTl, 



M6 



9M6 



a. 



(63) 



The four-quark structures appearing in the angle brackets, as well as the factor 
in front of them, are normalized at Q. The anomalous (logarithmic) dimension of 
this particular combination nearly vanishes [|I|; we will ignore it in evolving the 
four-quark operator at hand down to fi. Applying factorization, as was explained in 
Sect. 3, we get for the four-quark term in /npc(M^) 



4487r3 
'81M6 



tts(^)(gg(/i))^ 



M6 



0.03x4. GeV^ 



(64) 



where, as previously, we have introduced a factor X4q to allow for possible theoretical 
uncertainties. Under the standard choice of parameters this factor is equal to 
unity. If we are on the right track it is reasonable to expect that x^g is close to unity. 

A huge amount of work was carried out in calculating two-loop corrections in the 
coefficients of various condensates. The results can be inferred, say, from Ref. HH 
or from the reviews EH] . Unfortunately, the existing uncertainty in the condensates 



themselves precludes us from taking advantage of the precision achieved in pertur- 
bative corrections. Therefore, we will limit ourselves to the leading order results for 
Cg and C4q. 
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Now, the stage is set and we can finally present the sum rule for I{M'^), 



4 /M^ 

I{M^) = 1 + - a I — 1 {1 + 0.729 



9 



e"^ 





- 0.386 





















+ 



0.14X4g 



(—)' 



(65) 



where M in the second line of Eq. (|65|) is measured in GeV. The term M""^ is 
due to the gluon condensate, while the term is due to the quark condensate. 
With the "standard" numerical values of the gluon and quark condensates in the 
the nonperturbative part xq = x^q = 1; the dimensionless constants xg and x^q 
allow the gluon and four-quark condensates to "breathe" to a certain degree. For 
instance, > 1 would imply that the actual value of the gluon condensate is larger 
than the "standard" value, and so on. 

We are ready to examine the situation in the p-meson channel. Figure 16 shows 
experimental data on the spectral density in the vector isovector channel, measured 
in the r decays. The data points are from Ref. |^^. Note a remarkably close 
resemblance with our toy model number two. If we take the beginning of Fig. 8 
and expand it to make the energy scales coincide the curves will essentially repeat 
each other. There exist a few extra points in the r decays, spanning the interval 
of s from 2.7 to 3 GeV^, and some data above 3 GeV^ from the e~^e~ annihilation, 
but the corresponding error bars are so large that plotting these points would just 
obscure the picture. 

There is a subtle point which I have to discuss here. Directly measurable in the 
hadronic r decays is the sum of the vector and axial spectral densities. To obtain 
the spectral densities separately one has to sort out all decays by assigning specific 
quantum numbers to each given final hadronic state. In the majority of cases such 
an assignment is unambiguous. For instance, two pions (whose contribution is the 
largest) can be produced only by the vector current. Some processes, however, can 
occur in both channels, for instance, the KKtt production. Using certain theoretical 



arguments it was decided in Ref. that around 3/4 of all KKir yield should be 
ascribed to the vector channel. Other theoretical arguments [Q, which seem to 
be much more convincing to me, tell that virtually all KKtt production should 
take place in the axial channel. Therefore, in dealing with the data, I will subtract 
the KKir yield from the ALEPH data points. I hasten to add that the subtracted 
quantity is small (see Fig. 16), and this subtraction is not very essential for a general 
picture I draw here by broad touches. 

The solid curve in Fig. 16 is a best fit by a smooth curve representing the sum of 
two Breit-Wigner peaks (the first one is a modified Breit-Wigner taking into account 
threshold effects important for the p meson). The dashed curve is what I believe 



■'Warning: the energy scale (the horizontal axis) in Fig. 8 presents s in the units of while 
that in Fig. 13 is in GeV^. 
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0.5 1 1.5 2 2.5 



s, GeV ^ 

Figure 16: The spectral density in the vector isovector channel measured in r decays. 
The data points belong to ALEPH. The solid curve is a best fit. The dashed curve 
is the best fit minus KKtt (the latter contribution is shown by the dotted curve) 
which, I think, should not have been included in the vector spectral density, see 
text. 
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p (s) versus s, in GeV ^ 



perturbative asymptotics 



1 2 3 4 5 6 7 

Figure 17: "Theoretical experimental" spectral density in the p meson channel. The 
solid curve presents experiment (see text), while the dashed curve is the perturbative 
result for /3(s), including terms up to 0(af), with A = 0.2 GeV. 

the actual spectral density is. The data points above 1 GeV^ carry a noticeable 
uncertainty. Since we are going to use the spectral density only in the integrals, 
where many data points are summed over, these individual errors can be neglected, 
since they will be statistically insignificant in the integrals. Systematic uncertainty 
might be important, but since nobody knows how to estimate it, I will ignore it for 
the time being. For our limited purposes we can consider the dashed curve in Fig. 
16 as an exact experimental result for the spectral density. 

How does the experimental spectral density might look above 2.7 GeV^? This 
question is irrelevant for calculating the p-meson parameters. It is very instructive, 
however, to have a broader perspective. I will show you what I call a "theoretical 
experimental" spectral density, and then will make a short digression explaining 
where the answer comes from. 

My expectations are summarized in Fig. 17. The curve to the left of the arrow 
is just the fit to the experimental data, as explained above (it identically coincides 
with the dashed curve in Fig. 16.) To the right of the arrow I have attached a 
"tail" which approaches a smooth asymptotic prediction for p(s) (the dashed line) 
given in Eq. (|59D, in an oscillating manner. The tail and the experimental data are 
smoothly matched at 2.65 GeV^. 
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Why I think that the actual spectral density, when and if it is measured, will 
approach the smooth asymptotic curve with oscillations rather than monotonously? 
The reason is of a very general nature |TB|. In a nut shell, there are exponential 



terms that are not seen in the truncated OPE series for 11 (Q^) in the Euclidean 
domain. When these terms are analytically continued to the Minkowski domain, 
where Im n(s) belong, they show up as oscillating terms. The smooth asymptotic 
prediction is obtained from the truncated series and, thus, carries no traces of the 
exponential/oscillating terms. We will return to the issue later on (see Sect. 5). 

Although the very existence of oscillations in the spectral densities may be consid- 
ered firmly established, the question of how rapidly they die off is highly non-trivial 
and has no unambiguous answer in today's theory. In resonance-inspired models the 
damping factor is ~ exp(— const s) (see Sect. 2.2), in the instanton-based models 
Il48| it is ~ exp(— const a/s). The tail added in Fig. 17 is taken to be proportional to 
exp(— TTi/i) cos(7rs + const). This rather eclectic formula is chosen for sophisticated 
reasons which need not concern us here. The task which we address is illustrative, 
anyway. I certainly cannot guarantee that around 3.5 GeV^ the value of p{s) is 0.97, 
as it is shown in Fig. 17, but I could bet that a shortage of the spectral density in 
this domain will be observed in precision measurements, so that p{s) at 3.5 GeV^ 
will lie somewhere around unity. I emphasize again that the precise form of the tail 
does not affect the sum rule calculation of the p meson parameters, and is discussed 
only for completeness of the picture. 

With the experimental spectral density in hands, we can calculate the "experi- 
mental" value of the integral /(M^). The corresponding result is presented in Fig. 
18 by the solid curve. 

The agreement is excellent; it is even better than one could expect apriori. We 
see that the conspicuous p peak and all further twiddles characteristic to p{s) are 
washed out. If we descend from larger to smaller values of M, the behavior of /(M) 
is flat down to M ~ 0.8 GeV, i.e. down to the p meson mass. At M ~ 0.8 GeV 
the regime smoothly changes, the curve dives down, and at still lower values of M, 
approaches the exponential asymptotics (this domain is not shown in the figure). 
The boundary value of M where the regime changes is correlated with the mass 
of the lowest state, the p meson in the case at hand. This is practically the only 
characteristic dimensional parameter in I{M). 

Let us pretend now that we do not know the spectral density and want to use 
the sum rules to determine the p meson mass. It is quite clear that the truncated 
condensate series does not allow one to go to the limit M — > where this determi- 
nation would be exact: the series explodes. We must stay inside the window where 
the expansion is still under control. Correspondingly, our determination of Mp is 
going to be approximate. The strategy can be formulated as follows. Let us assume 
that somebody shows us a sketch of the spectral density presented in Fig. 10 with 
the numbers along the horizontal and vertical axes erased. This sketch is known to 
correctly reproduce the basic features of the actual spectral density but leaves us 
ignorant as to where the p-meson peak lies and what is its height. We insert this 
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p meson sum rule 



Figure 18: I versus (in GeV^) in the p meson channel: confronting experiment 
and theory. The theoretical curve corresponds to xq = 0.8, x^q = 1.3 and A = 0.2 
GeV. 
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sketch in our sum rule, do the integral 



7rM2 



ds p{s)e 



fit the numbers in such a way as to be as close to the theoretical prediction (|65[) 
inside the window - click, click - the scales are restored, and there come out the 
p mass and the coupling constant. Since inside the window the p meson saturates 
the integral at the level ~ 90% , the fact that our continuum model is a caricature 
(minor details are lost) is unimportant. Even if we are off by a factor of two in this 
model this will affect our estimates referring to the p meson at the level of ~ 10%. 

This example is quite typical. A similar situation takes place for all classical 
low-lying hadrons: those built from the light quarks, heavy, and light and heavy. 
I do not define here precisely what the "classical meson" means but will return to 
this point later, after considering some nonclassical channels. I must admit that the 
p meson is the example where the SVZ method demonstrates its best facets. The 
window is sufficiently broad, everything is clean. As we proceed to higher spins, 
the mesons become larger in size. A snapshot of a high-spin meson would show a 
string-like picture, with a longitudinal size of the "sausage" much larger than its 
transverse size. Under the circumstances one could hardly expect that the SVZ 
method would work. And, indeed, it ceases to be informative for spins higher than 
two El. 



5 Basic Theoretical Instrument — Wilson's OPE 

The theoretical basis of any calculation within the SVZ method is the operator prod- 
uct expansion (OPE). It allows one to consistently define and build the (truncated) 
condensate series for any amplitude of interest in the Euclidean domain. The phys- 
ical picture lying behind OPE was described above: consistent separation of short 
and large distance contributions. The former are then represented by the vacuum 
condensates while the latter are accounted for in the coefficient functions. Thus, the 
operator product expansion in the form engineered by Wilson is nothing but 



a book-keeping procedure. Wilson's idea was adapted to the QCD environment in 
Ref. PI. 



Although OPE is used in an innumerable amount of works since the mid-1970's, 
there are no good text-books or reviews devoted to this issue. Moreover, one typi- 
cally encounters a lot of confusion in the literature. Below I will try to explain the 
Wilsonean procedure using, as an example, the T product of the two vector currents 
defined in Eq. (0). 

Technically it is more convenient to deal with the D function (^) rather than 

D{q') = Y.Cn{q;p){OM) (66) 
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Figure 19: A typical multiloop graph contributing to the D function. The number 
of such graphs grows factorially with the number of loops. 



where the normalization point /i is indicated explicitly. The sum in Eq. (66) runs 
over all possible Lorentz and gauge invariant local operators built from the gluon 
and quark fields. The operator of the lowest (zero) dimension is the unit operator I, 
followed by the gluon condensate G^^,, of dimension four. The four-quark condensate 
gives an example of dimension-six operators. 

At short distances QCD is well approximated by perturbation theory. The co- 
efficient functions Cn absorb the short-distance contributions. Therefore, as a first 
approximation, it is reasonable to calculate them perturbatively, in the form of ex- 
pansion in as{Q) ~ (InQ)^^. This certainly does not mean that the coefficients Cn 
are free from non-perturbative, non-logarithmic terms of the type ~ Q~'^ where 7 
is a positive number. 

To substantiate the point let us consider the coefficient of the unit operator. 
The Feynman graphs for D{q^) of the lowest order are depicted in Figs. 5 and 
13. Assume that the momentum flowing through the graph is Euclidean and large, 
Q 00. If all subgraphs in these and similar graphs are renormalized at Q, the final 
result, being expressed in terms of the running coupling constant dsiQ), is finite. 
The virtual momenta saturating the loop integrals scale with Q as the first power 
of Q. At any finite order the perturbative series is well-defined. 

At the same time, if the number of loops n becomes very large, in a random 
graph (Fig. 19) each loop carries momentum of the same order of magnitude, and 
the total momentum Q is shared between many lines, so that the characteristic 
virtual momentum is proportional to Q/n. The contribution of each graph to Cj is 
of order [^^(Q)]'^, but the number of distinct diagrams u grows factorially, u ~ n\, 
so that the expansion in as{Q) is asymptotic (for an exhaustive discussion of this 
phenomenon see Ref. jSTH). 

Note that this factorial divergence has nothing to do with the renormalon di- 
vergence considered in Sect. 3. There, we had one specific graph of the n-th order 
whose contribution was proportional to n\. This factorial was a spurious artifact 
of an improper calculation where the domain A; ~ A was included. Introducing a 
lower cut-off at yU = several units times A and discarding the contribution coming 
from fc^ < /i^, as it should be done in the consistent calculation of Ci{fi), we would 
eliminate the renormalon divergence. At the same time, the factor n\ counting the 
number of graphs of the n-th order certainly cannot be eliminated in this way. 

If n is not too large, so that Q/n > /x, the graphs of Fig. 19 present a legit- 
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imate contribution to Ci{^). It is quite clear that the tail of the series with 
n > l/as{Q) ~ InQ generates a non-perturbative contribution of the type 




AC I ~ exp {-2n^/asiQ)) ~ TT • (67) 



Here 7 and n are numerical constants. An explicit example of such terms is provided 
by the so called direct instantons ||^. Integration over the sizes p of the direct 
instantons is saturated at p ~ Q^^. The constant 7 = 1 for one instanton, 2 for two 
and so on. 

The correspondence between the Feynman graphs with l/as{Q) loops and the 
small-size instantons is not straightforward. I will not go further in this issue, 



referring the interested reader to Ref. fST]]. At the qualitative level the inevitability 



of occurrence of the hard non-perturbative terms ([67| ) is evident. At the moment 
the only semi- quantitative framework we have at our disposal for their evaluation is 
the instanton mechanism. 

Equation (^) gives a non-condensate power term. Generically we will call such 
terms hard non-perturbative. They may or may not be important numerically de- 
pending on the particular value of Q under consideration. The value of 67 need not 
be integer, generally speaking; numerically it is very large. This means that once 
we cross the boundary where the argument in the brackets in Eq. (^) becomes less 
than one, Q > nA, these terms immediately become totally unimportant. On the 
other hand, below this boundary they are so large that no expansion is possible. 

It is assumed that inside the window the terms (^) are unimportant; they are 
completely disregarded in the practical version of OPE (Sect. 6) which, thus, applies 
to the values of Q above the critical point. Note that the large value of 67 results 
in a large degree of suppression of such terms after borelization. 

Although the assumption above works very well in many "classical" channels, 
it seems to fail in exceptional cases of "non-classical" mesons. This issue will be 
addressed in more detail in Sect. 8. 

I pause here to make a remark concerning non-dynamical power terms in C/(/i) 
whose origin is related to the introduction of the normalization point p. Let us 
return to the one- loop graph of Fig. 13 and the corresponding expression (^21). One 
should not forget that in doing the loop integrations in Cj{fi) we must discard the 
domain of virtual momenta below n, by definition. Subtracting this domain from 
the perturbative loop integrals we introduce power corrections of the type (/i^/Q^)" 
in C/(/i), by hand. For instance, in the first order in as 



as{Q) as p 



2 



+ ^W/x)) + ..., (68) 



where I combined Eqs. (E3) and (pTf). The expression in the square brackets is the a 



part of the coefficient Ci{fi). The explicit /i dependence of this expression conspires 
with an implicit /i dependence residing in {Oc{p)) to ensure the p independence of 
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the physical quantity (i.e. the D function). If higher-order terms in the k"^ expansion 
of the F function (^31) were taken into account, we would get higher powers of /i/Q in 
the square brackets, which then should have been combined with higher condensates, 
e.g. GD^G and so on. 

If /i can be chosen sufficiently low, the explicit and implicit /i-dependent terms 
in Eq. ( |68D may be numerically insignificant and can be ignored. 

The coefficient functions in front of other operators generally speaking have the 
same structure as C/. The general situation is quite similar. For higher dimension 
operators it may happen (and, in fact, happens) that some of the loop integrations 
are not saturated at Q or /i but are, rather, logarithmic. Logarithms [ln((5//i)]''' 
occurring in this way are associated with the anomalous dimensions of the operators 
at hand. 

Let us have a closer look at the graph depicted in Fig. 13. As we already know, 
some of the lines can be in a special regime - the corresponding virtual momenta 
are soft and their off-shellness does not scale with Q. If the gluon line is soft, this 
gives rise to the gluon condensate, if some of the quark lines are soft we deal with 
the four-quark condensate. Graphically one distinguishes the soft lines by cutting 
them; then the remainder of the graph shrinks to a point while those lines that are 
cut form a local operator. The question is: what happens if we pass to the multiloop 
graph of Fig. 19 and start cutting more and more lines? 

Naturally, in the power series we proceed to the operators of higher and higher 
dimension. When the number of the soft lines becomes very large (of order of Q/A) 
it is conceivable that there are no hard lines at all: the external momentum Q is 
transferred from the initial to the final vertex through a very large number of quanta 
(growing like a power of Q), and none of the quanta carries momentum scaling with 
Q. Of course, in this situation one can not speak of individual quanta, one should 
rather use the language of typical field fiuctuations transmitting the momentum Q 
without having any Fourier components with frequencies of order Q. Contributions 
of this type are not seen in OPE truncated at any finite order. They correspond to 
the high-order tail of the condensate series which, analogously to the as series, is 



factorially divergent. In this way we obtain the so-called exponential terms [|T^ . 

The fact that the condensate series is factorially divergent in high orders is rather 
obvious from the analytic structure of the polarization operator n((5^). In a nut 
shell, since the cut in 11 (Q^) runs all the way to infinity along the positive real 
semi-axis of g^, the 1/Q^ expansion cannot be convergent. The actual argument 
is more subtle than that, but I would not like to go into details here referring the 
reader to Ref. where a careful consideration of the issue is carried out. The final 
conclusion is perfectly transparent. It is intuitively clear that the high-order tail of 
the (divergent) power series gives rise to exponentially small corrections ~ exp(— Q*^) 
where a is some critical index. 

The numerical value of a is correlated with the rate of divergence of high or- 
ders in the power series. At the moment very little is known about this rate from 
first principles, if at all. The best we can do is to rely on toy models. The sim- 
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plest example is again provided by instantons. This time one has to fix the size 
of the instanton p by hand, p = Pq. Then the fixed-size instanton contribution is 
C(exp(— Qpo))- The exponential factor is the price we pay for transmitting the large 
momentum Q through a soft field configuration whose characteristic frequencies are 
or order pg^. 

One can look at this example from a broader perspective. Consider the polariza- 
tion operator 11 in the coordinate rather than the momentum space. If the regular 
terms in OPE are in one-to-one correspondence with singularities of n(x) at a; = 0, 
the exponential terms, invisible in the truncated OPE, are related to the singulari- 
ties of Il{x) located at a finite distance from the origin [Q; the distance from the 
origin plays the same role as the instanton radius Pq. For our illustrative purposes it 
is sufficient to consider the singularities of the simplest possible structure, namely, 

ln(x^ + PQ), (x^ + po) In (x^ + Pq), and so on. (69) 



x'^ + pI' 

The Fourier transforms of these expressions have the generic form 

(gpo)-"i^n(gpo), n = l,2,... (70) 

where Kn is the McDonald function. At large Euclidean the corresponding 
contribution dies off exponentially. 

An = (Qpor'Kr^iQpo) oc (gpo)-"-'/' e-«^« , (71) 

in full accord with intuition regarding transmitting a large momentum through a 
soft field fiuct nation. 

The terms that are exponentially small for Euclidean values of become oscil- 
lating upon analytic continuation to the Minkowski domain, Q iQ. Analytically 
continuing AH given in Eq. ( |7TD we arrive at 

ImAn = (-l)"+i|(^po)-"J„,(v/^po) cx (-l)"+i^(V^po)— V2eos(v^po - 5„) , 

(72) 

where Jn is the Bessel function. The oscillating character at > is the most 
important signature of the exponential terms representing the high-order tail of the 
power series. 

Having fixed po is unrealistic. We would do much better than that by allowing 
Po to vary around a typical value of order of A~^. In other words, any reasonable 
model will introduce a distribution over po, with a weight function w{po) peaked at 

(ImAn) oc / {^spo^'UV^PoMpo)^ (73) 
J Po 

where the angle brackets denote po smearing. The weight function w{po) is model- 
dependent. Particular details of the result will certainly depend on w(po), but 
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the general pattern - the exponential fall-off of ImAII modulated by oscillations - 
remains the same under any reasonable choice of w{pq). This is the reason why the 
experimental spectral density in Sect. 4 is extrapolated beyond the range of the 
direct measurements in an oscillating mode (Fig. 17). The exponential suppression 
factor at large E takes the form 

expi-kE"") , (74) 

where the critical index a in various models lies in the interval < o" < 2. 

It remains to be added that the exponential terms are usually ignored in the sum 
rule analyses one encounters with in practical applications. 

6 Practical Version of OPE 

In practical calculations one can hardly go beyond several lowest-order terms in 
the as expansion of the coefficient functions and in the condensate expansion. For 
instance, in the p meson channel, which is the most advanced, four terms of the 
perturbative series are known in the unit operator. 

Thus, the approximation we make at the very first stage of building OPE is 
truncating both the perturbative and condensate series. By doing so we expect that 
those few terms that are kept represent -D(Q^) or J(M^) inside the window with 
sufficient accuracy. 

Besides truncating the series, one usually ignores the hard non-perturbative 
terms in the coefficient functions. The reason why these terms are neglected is 
quite obvious: no framework allowing one to reliably calculate them was worked 
out so far. From the instanton studies we know that such terms must have an 
extremely steep M dependence, so that above a critical value of M their impact is 
expected to be totally negligible. It is assumed that the critical value of M, where 
the hard non-perturbative terms die off, lies to the left of the window. 

This is not the end of the story, however. The third key element inherent to 
the practical OPE is a simplified treatment of the p dependence in the coefficient 
functions and the condensates. The perturbative parts of the coefficient functions 
are often borrowed from earlier calculations of relevant Feynman graphs which, 
sometimes, date back to QED. These calculations make no distinction between small 
and large virtual momenta. Moreover, in analytic QCD it is technically not always 
easy to explicitly carry out the Wilsonean separation of virtual momenta (higher 
than fi and lower than /i). At the very least, such a separation requires dedicated 
analyses. 

This is an obvious task for the future. A simplified strategy of today is as follows. 
Assume that we do not plan to go beyond the one-loop correction of Fig. 13 in the 
perturbative part of Cj and beyond the leading-order result for Cq (and higher 
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condensates). The corresponding proper expression for the D function is 



D=l + aJ deF{e,Q') + ^—{OG{fi)) (75) 

(cf. Eq. (HD). Here F' is dF/dP at P = 0. The dimension -six condensate GD^G 
enters with the coefficient proportional to F'' at k"^ = 0, and so on. 
Now, let us add and subtract to the right-hand side the integral 



a. 



r dk^Fik^Q""). (76) 
Jo 



This integral looks exactly as the one-gluon contribution naively continued below 
/i^. Of course, below /i^ the gluon propagator has nothing to do with l/k"^, but since 
the transformation is identical, we do nothing wrong. Equation (|75| ) then takes the 
form 



D = l + aJ deF{k\Q^) + ^— (Og(/x))-^ / k'^dk 

Jo 6 TT-^ Jo 



,2 7,2 



(77) 



where in the square brackets, I expanded the integrand in k"^ and kept only the 
leading term. The next term proportional to k"^ would be important at the level of 
dimension-six condensate, etc. 

The first integral on the right-hand side is the full perturbative one gluon cor- 
rection to Cj, with no Wilsonean separation (it is equal to as/ir). The term in the 
square brackets can be called an effective one-loop gluon condensate. It is obtained 
from the genuine condensate by subtracting from it its perturbative one-loop expres- 
sion. By construction, the fi dependence of the effective one-loop gluon condensate 
cancels provided we do not go beyond the accuracy specified above. 

Adding the integral (^) to the right-hand side we achieve the desired goal. Now 
the coefficient Cj, to order 0{as), is given by full perturbative graphs of Fig. 13. 
To avoid double counting we subtract the same contribution from the condensates. 
More exactly, we subtract almost the same contribution. Since the sum over con- 
densates is truncated at some finite order uq, the best we can do is to subtract from 
the condensate part the integral 

no 1 

as / de y ^F(")(fc2 = 0, Q2) (7g) 
Jo ^1 n\ 

In this way we subtract from each condensate its "one-loop perturbative value" . 

It is quite clear that this strategy of converting the genuine Wilson OPE into the 
practical version is admittedly approximate and works only if we limit ourselves to a 
given number of loops and given number of terms in the condensate expansion. The 
effective condensates obtained in this way should be supplemented by a superscript 
indicating the number of loops in the subtracted part. Their numerical value is fi 
independent but does depend on the number of loops. 
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This approximate procedure (i) leads to a loss of factorization of short and large 
distance contributions inherent to Wilson's OPE; (ii) can not be systematically 
generalized to arbitrary number of loops since there is no way to unambiguously 
define the integrand to all orders in the small fc^ domain (iii) only approximately 
avoids double counting due to the necessity of truncating the condensate series at 
a finite order (the lower the order we truncate the larger the error). Moreover, 
the effective condensates are not universal, strictly speaking; they become process- 
dependent. If this non-universality were numerically significant, the SVZ method, 
as it exists now, would be undermined. 

Therefore, the practical version is useful in applications only provided /i^ can 
be made small enough to ensure that the "one-loop perturbative" contributions to 
the condensates are much smaller than their genuine values and, at the same time, 
as(/i^)/7r is small enough for the expansion to make sense. The existence of such 
"/i^ window" is not granted a priori and is a very fortunate feature of QCD. We 
do observe this feature empirically although the roots of the phenomenon remain 
unclear. 

Although the origin of the /i^ window, ensuring the applicability of practical 
OPE is not understood, one can find similar situations in simplified settings. The 
best example of which I am aware is the two-dimensional 0{N) sigma model in the 
limit of large N. In this limit, the model is exactly solvable; therefore any question 
of interest can be exhaustively answered. The model bears remarkable parallels 
to QCD - it is asymptotically free, and the mass scale is generated dynamically. 
A surprising result was established some time ago |Q. In the large N limit all 



vacuum condensates in this model are fi independent, and the practical version of 
OPE becomes exact; fi dependence in the condensates and in the coefficient functions 
comes only at the level of 1/iV. If iV is large these terms are parametrically small. 
(It would be great to find a hidden parameter which would explain the suppression 
of the fi'^/Q'^ terms in QCD in the same way 1/N explains it in the sigma model.) 

It is often asked whether it is legitimate to retain in the truncated expansions 
logarithmic and power terms simultaneously Indeed, formally any 1 / (In Q)" term 
is parametrically larger that, say, l/Q^ (moreover, any power term is parametrically 
larger than the exponential ones). If so, the theoretical uncertainty due to truncation 
of the logarithmic series is formally larger than even the lowest-order condensate 
correction. 

The logarithmic and power terms in practical OPE have distinct physical nature. 
In some instances power terms describe effects that do not show up in perturbation 
theory. This is valid in all cases where the chiral quark condensate is involved. In 

^^Within practical OPE we automatically forbid to ourselves questions concerning any aspects 
of the high order behavior. 

^^A natural abbreviation for practical OPE would be POPE, but I do not risk to put this 
abbreviation into circulation. 

^•^In discussing the issue of duality violations one deals with the logarithmic, power, and expo- 
nential terms simultaneously, see Ref. Es] . 
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other cases, when the power terms are not different in their structure from those 
occurring in perturbation theory, there is a strong numeric enhancement of the 
power terms (see Sect. 3). This mysterious fact - the numerical enhancement of the 
vacuum condensates - explains why Wilson's OPE can be substituted in QCD by 
the practical version, at least in the classical channels . 

In the early days of QCD, in 1970's, it was common to assume that the coefficient 
functions in the operator product expansion are saturated by perturbation theory, 
with no separation of virtual momenta. Although the difference between Wilson's 



OPE and the practical version was realized long ago ||T8[, very little has been done 
to investigate this difference and perfect the procedure. Only now we are witnessing 
attempts in this direction, mostly in connection with the heavy quark theory (e.g. 

1,00)- 



7 Low Energy Theorems 

Although the low-energy theorems are independent of the sum rules, at least some 
of them were derived in connection with the development of the SVZ method, and 
were later combined with the sum rules giving rise to an intriguing observation of 
the hadronic non- universality. This observation is the topic of the next section. Here 
we will dwell on the low-energy theorems related to the trace anomaly of QCD. 

In the limit of massless quarks (and we will never leave the limit = 0) 
the classical Lagrangian of QCD is free from dimensional parameters. The scale 
parameter appears at the quantum level, 

where Mq is the ultraviolet cut off, go is the bare coupling constant, and I ignored the 
second and higher-order terms in the f3 function. (The corresponding modifications 
are quite trivial and are suggested to the reader as an exercise.) 

The occurrence of the scale parameter (^) is in one-to-one correspondence with 
the trace anomaly (|39D . Let us denote the trace operator 

e^^^a. (80) 

Then, it is not difficult to show that for arbitrary local operator O 

lim |z J d^a;e^^^'(T{O(0),a(x))c| = -d (O) , (81) 

where d is the (normal) dimension of the operator O, and the subscript c denotes 
the connected part (it will be omitted hereafter). If 

C = a oc G'^ 
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then the dimension c? = 4, so that 



(82) 

For the quark operator qq the dimension (i = 3, and so on. 

The derivation is straightforward. First, we rescale the gluon field, = G^^, 
so as to factor out the go dependence in the Lagrangian. Instead of Eq. (|) we now 
have ^ 

£ = --^ Q%Q% + quark term . 

Then we observe that 

^ j d'x{T{0{Q\G\x)]) = -^(Y^ {O). 

Since {O) = A'^, we use Eq. (^) to perform differentiation on the right-hand side. 
In this way we arrive at the formula 

In fact, the simple derivation outlined above should be supplemented by some 
regularization since the matrix elements (O), as well as the two-point functions (plj), 
are divergent unless the cut-off at /i is introduced. In Ref. (see Appendix B) 
it is shown that the regularization procedure does not affect Eq. ( |HT| ) provided the 
cut-off is introduced in a concerted way. 

Denote the two-point function on the left-hand side of Eq. ( P^ ) as Us^Q"^); the 
subscript s marks the scalar glueball channel. If ns(0) is known, this knowledge can 
be immediately traded for a power term in the Borel-transformed sum rule. Indeed, 
in this case, instead of borelizing the dispersion relation for Ils{Q^), one can deal 
with the dispersion relation for 

ns(g^) - uM 

After borelization one gets a prediction for the integral 

1 r ds 



7rM2 



/ -Imn,(.)e-^/^' 
J s 



Note the occurrence of an extra factor s ^ in the weight function in the integrand 
(cf. Eq. (^). At the very end 

n.(o) = ^(Oc) 

is transferred to the right-hand side of the sum rule where it acts as a coefficient in 
the corresponding power correction. 
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In this way one arrives at 



1 + 



IGtt^ 1 .^^ , 



hal{M) M4 



(83) 



where the first term in the square brackets corresponds to the free gluon loop, 
the second term is due to ns(0), while the ellipses stand for all other corrections, 
perturbative and non-perturbative. Comparing the gluon condensate correction to 
that in the p meson sum rule, Eq. (^), we observe, with astonishment, a very 
strong numerical enhancement of the relative coefficient, by a factor 
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The standard condensate expansion, with a few first terms kept, does not match 
this huge low-energy constant. Attempts to smoothly match this gigantic power 
term with other non-perturbative corrections and to saturate the sum rule in the 
scalar glueball channel led us [0] to the conclusion of extremely violent distortions 
of the spectral density up to the scale Sq ~ 25 GeV^, to be compared with sq ~ 1.5 
GeV^ in the p-meson sum rule, see Eq. ( PUD and Fig. 10. The scale up to which 
strong violations of asymptotic freedom extend in the scalar glueball channel is 
unconventionally large. 

Thus, a new (numerically) large scale was discovered, that does not show up in 
the classical channels. It was predicted to have a major impact on the properties 
of the "non-classical" mesons. The physical picture lying behind this phenomenon 
is discussed in Sect. 8. The quantum numbers of the mesons affected by the "su- 
perstrong" interaction with the vacuum medium were found to be in one-to-one 
correspondence with the presence of the direct instantons. This observation served 
as an initial impetus for the development of the instanton liquid model of the QCD 
vacuum ||5^ (for a review see ||56[|). 

The topic of the QCD low-energy theorems is very vast; it reaches out to such 
profound issues as the subtleties of the 6 dependence, the rj' problem, effective dilaton 
Lagrangians, the 't Hooft matching condition, etc. It has never been reviewed in 
full, although it certainly deserves attention. This endeavor goes far beyond the 
scope of the present lecture, and I will limit myself to fragmentary remarks on the 
literature. 

The low-energy limit of the n-point functions generated by GG was studied in 
A celebrated mass formula for rj' was obtained in An elegant low- 

energy theorem for the coupling of the Goldstone mesons to the gluonic sources 
was derived in Ref. [^]. The couplings of photons to the gluonic sources were 
elaborated in Refs. |61|. Low-energy theorems instrumental in the sum-rule 
analyses, additional to those considered in P], were found in [|^. Selected aspects 
of the subject are reviewed in p3|. 
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8 Are All Hadrons Alike? 



The first, superficial, impression of the hadronic family may leave one quite bored. 
There are no obvious small or large parameters, all quantities are of order one in the 
appropriate units, and nothing special attracts attention. In particular, one might 
think that the characteristic hadronic scale is given in all cases by the slope of the 
Regge trajectory, i.e. universal a' ~ IGeV^^. The masses of the "most" typical 
hadrons, nucleons and p mesons, are of this order of magnitude. Roughly at the 
same momentum transfers scaling sets in in deep inelastic scattering heralding the 
onset of asymptotic freedom - the observer sees (almost) free quarks with modest 
corrections due to the gluon exchanges. 

A careful reflection shows, however, that the hadronic family is far from being 
that monotonous. Messages of its non-universality have been coming from various 
sides, although it was not so easy to decipher them. 



Many years ago Witten posed a question ||58[| : 

"There is an obvious troublesome question. If m^, ^ l/N^ and the l/N^ expan- 
sion is usually a good approximation, why is m^/ not much smaller than its actual 
value of almost one GeV? 

This is not merely a problem oil/Nc expansion. It is a more general phenomeno- 
logical problem. The rj' — ir mass splitting violates Zweig's rule, since it involves qq 
annihilation. Since Zweig's rule is usually rather good, why is the rj' — it splitting 
so large?" 

A number of the glueball candidates have been experimentally observed in the 
last decade. Still, none was unambiguously proved to be a glueball. Let us ask 
ourselves: 

Why, unlike the quark mesons, the masses of the glueballs are unusually heavy 
or they have an unusually large mixing with the quark states which makes their 
widths large enough to escape clear experimental identification, in defiance with the 
1/Nc counting? 

The answers to all these questions naturally come with the observation ||^ of large 
scales in the non-classical channels. Yes, m^, and Zweig's rule violating parameters 
are suppressed by l/N^ and are small, but they are small compared to the intrinsic 
scale So ~ 25 GeV^, rather than 1 GeV^. 

The existence of the diverse scales in the hadronic family is due to the fact that 
various currents injecting the valence quarks/gluons in the vacuum medium interact 
with this medium non-universally. If in the p meson (and similar classical mesons) 
the impact of the vacuum fields on the valence quark-antiquark pair is modest and 
is reasonably well described by the average densities of the vacuum fields, in the 
glueballs the coupling to the vacuum fields is much stronger, so that they actually 
feel fine details (a grain structure) of the vacuum medium. The strongest coupling 
to the vacuum fields takes place in the channels with the total spin zero, both quark 
and gluon (scalar and pseudoscalar quarkonia and glueballs). Phenomenologically 
deviations from Zweig's rule are most drastic here. 
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The occurrence of the "superstrong" couphng manifests itself in an abnormally 
large critical scale which may or may not materialize as the mass of the lowest-lying 
state with the given quantum numbers. By critical scale I mean the energy at which 
the asymptotic perturbative regime sets in in the spectral density. 

A combined analysis of the sum rules and the low-energy theorems has led us [Q] 
to a number of qualitative (and, in many instances, semi-quantitative) predictions. 
First of all, the hierarchy of glueball masses was established. The scalar glueball was 
predicted to be the lightest, with mass M(0"*") ~ 1.5 GeV (see |@]). Its coupling to 
the two-pion state was shown to be very large, with no trace of the Zweig suppression. 
The pseudoscalar and tensor 2+ glueballs were predicted to lie in the vicinity of 2 
GeV. 

In both scalar and pseudoscalar channels the valence gluons injected in the vac- 
uum are affected by the vacuum medium in the most drastic way, and the critical 
scale is the largest. In the tensor glue channel the critical scale was expected to 
be smaller than in the scalar/pseudoscalar channels, although larger than for the 
classical mesons. 

Why the 0^ channels are so special, with the strongest reaction of the medium? 
It was natural to assume that the phenomenon had to do with the structure of the 
dominant vacuum fluctuations. It was noted |^ that instantons would do the job, 
since the direct instantons show up exactly in the 0"*^ channels. This was sufficient 
for a qualitative classification. Quantitative instanton-based models of the QCD 
vacuum appeared later ||5^, they automatically incorporate all relevant features 
of this peculiar picture. Hunting for instantons in the vacuum is the vogue of the 
day in the lattice community. There are indications that they are abundant and 
almost saturate the string tension. 

It remains to be added that, many years after these predictions, we are witnessing 
how other approaches started recovering the very same pattern of non-universality. 
Recent calculations in the instanton liquid model produce - not surprisingly - a 
very close hierarchy of masses and mixing parameters [^]. What is interesting, this 
model explains where the large non-perturbative scale in the scalar glueball channel 
is hidden. The violent vacuum fiuctuations squeeze this state to an unusually small 
radius because of a very strong attraction. The radius of the 0"*" glueball in this 
model is only 0.2 fm, to be compared with 0.6 fm in the p meson case |^5[. The 



ratio of the radii squared is ~ 1/10. It is clear that no constituent model compatible 
with quantum mechanics can accommodate a state light and narrow simultaneously 
(see the end of this section). 

This part of the hadronic family is being intensively explored on the lattices 
too. Of course, the light quarks are usually non-dynamical on today's lattices, and 
a typical lattice site of 0.2 fm is comparable with the radius of the 0"*" glueball. 
With all these reservations in mind it is still instructive to compare what the lattice 
community learned about the non-classical hadrons. To this end I borrowed relevant 
numbers from the talks and original publications to be cited below: 

(i) The lightest glueball is scalar, and its mass is close to 1.6 GeV; 
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(ii) The tensor glueball is significantly heavier; its mass is shghtly above 2 GeV; 
the pseudoscalar glueball is found at approximately the same place, with larger 
errors; 

(iii) The sizes of scalar and tensor glueballs were found to be drastically different. 
This can be inferred from the different magnitude of the finite size effects (see 
e.g. 10), or seen directly, in the glueball Bethe-Salpeter amplitudes (or "wave 
functions") |^ The radius of the scalar glueball was found to be 0.2 fm, 
while that of the tensor one is much larger, approximately 0.8 fm. A similar lattice 
measurement of the p meson size yields approximately 0.5 fm. 

Needless to say that the lattice calculations provide us only with final numbers, 
with no insight as to the mechanisms of the non-universality. To the best of my 
knowledge, no lattice publication on this topic mentions the previous sum rule pre- 
dictions. True, the lattice results have a potential to be perfected in the future, 
while the accuracy of the sum rule method is admittedly limited. 

To make the highly nontrivial nature of our findings more contrast let us conclude 
this section by comparing them to the expectations of conventional models, say, the 
bag model. The quark and (electric) gluon modes in the spherical cavity have 
energies 2.04/i? and 2.7/i?, respectively, where R is the cavity radius. Thus, the 
glueball states are expected to be only marginally heavier than the quark states. 
In particular, if the spin-dependent forces are neglected the model predicts that 
M(2++) K. M(0++) 1 GeV and M(0-+) ^ 1.3 GeV. Including the spin-dependent 
forces changes these numbers rather insignificantly |j70|. The complex picture of the 



vacuum medium, which most strongly affects the scalar and pseudoscalar mesons, 
especially glueballs, is completely missed. 



9 Ecological Niche 

In the 1980's the original strategy of the SVZ method was tested, with a remarkable 
success, in analyzing practically every static property of all established low-lying 
hadronic states. The method was developed in various directions - three-point 
functions, inclusion of external electromagnetic and other auxiliary fields, light-cone 
modifications and so on - which allowed one to expand the range of applicability to 
such advanced problems as magnetic moments |7l|], form factors at intermediate 



momentum transfers, weak decays [|73[ and structure functions JT^ of deep inelastic 



scattering at intermediate x. At the initial stage the method was essentially unchal- 
lenged since the lattice calculations were lagging far behind burdened by multiple 
internal problem of lattice QCD. 

Now the situation has changed. The machinery of lattice QCD was bettered, and 
it became a powerful tool in many problems where qualitative insight is unimportant, 
and the prime emphasis is on numbers. Let us remember that the SVZ method is 
admittedly approximate: typically, theoretical accuracy is at the level of 20%. In 
some cases it does not work at all. 
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In order to survive, the approach based on OPE and the sum rules has to find a 
range of apphcations in the hadronic physics, where it could successfully compete. 
I must say that it copes with the challenge. Its strongest advantage is the analytic 
character of calculations which can be continued to the Minkowski space, term by 
term. Thus, the most promising directions of growth are those where the Minkowski 
kinematics is entangled in this or that way and/or the processes under consideration 
develop in more complicated media, rather than in the vacuum medium. (Such 
problems are very difficult for lattice QCD which operates with vacuum correlation 
functions numerically evaluated in the Euclidean space.) Let me give three examples. 

1) Particle decays with a large spatial momentum carried away by produced 
hadrons, for instance, B piu. If the invariant mass of the lepton pair is small, 

~ 0, the spatial momentum of the p meson is of order Mb/2 ^ A. Treating 
the problem in the Euclidean domain we will have to face the necessity of an ex- 
tremely remote analytic continuation. At the same time the light-cone sum rule, to 
be discussed in Sect. 10, are perfectly fit to deal with this problem. 

2) . Heavy flavor sum rules. They allow one to determine basic parameters of the 
heavy quark theory - the masses of the heavy quarks, form factors at zero recoil, 
and so on - in terms of directly measurable quantities. This is a close relative of the 
SVZ sum rules. The main distinction is that instead of the vacuum condensates one 
deals with the expectation values of various local operators over the heavy hadron 
state, for instance the B meson. This approach will be also considered in brief in 
Sect. 10. 

3) Pre- asymptotic effects in inclusive processes at high energies. A typical prob- 
lem from this class is determination of the lifetime differences in the b quark family. 
Although the task is somewhat different from those usually treated in the SVZ sum 
rules, the method of its solution is exactly the same. The key technical element is 
the operator product expansion. In the given context it is imperative to address, 
additionally, an issue going beyond OPE (more exactly, its practical version) - du- 
ality violations. The topic of the deviations from duality is of paramount practical 



importance; it came under renewed scrutiny recently |16 



In all these and other instances the OPE-based methods remain to be the only 
theoretical tool available on the market today. We see that the approach is suffi- 
ciently rich and fiexible to remain viable in the future. I am sure it will continue to 
play a key role in solving many applied problems in the hadronic physics in the next 
decade. The potential for further growth is there, the active stage is not yet over. 



10 New Developments 

To give you a fiavor of what is going on in the this field I will pick up two or three 
subjects which come to my mind first. My selection by no means presents the full 
picture. 
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The heavy quark mass l/rriQ expansions that flourished in the 1990's revived 
interest in the SVZ sum rules for hadrons containing one heavy quark |^ (in practice, 
b quark). It was realized - and this is a new element - that (i) the hybrid logarithms 
occur in the theoretical part of the sum rules, and (ii) they may result in a strong 
enhancement of the gluon radiative corrections. Hybrid logarithms depend on the 
ratio mq/fi, they were unknown before 1988. A typical example where summation of 
the hybrid logarithms significantly shifts the answer is Jb- I single out this problem 
because this parameter is fundamentally important in the heavy quark physics and 
because of an instructive story attached to it. Early sum rule calculations of Jb, 
which missed then undiscovered hybrid logarithms, yielded a number close to 130 



MeV ||Tg]. The main impact of the hybrid logarithms is a change of the argument 
of the running coupling constant in the sum rule: instead of as{mh) the radiative 
corrections, as it turns out |7^, are governed by as{l GeV). As a result, the sum 
rule prediction jumped up to /b = 160 MeV, with the theoretical uncertainty 20 to 
30 MeV 

Meanwhile, fs was the object of the multiple lattice studies. The first lattice 
measurements produced an unbelievably huge value, somewhere around 250 MeV. 
With elimination of various sources of the systematic uncertainties the lattice num- 
ber has been systematically decreasing. The modern number, quoted at recent 
conferences, is close to 180 MeV, with the error in the same ballpark as the sum rule 
uncertainty. Thus, the initial contradiction between the two methods faded away 
giving place to perfect agreement. 

Along with the technical developments of the type I have just mentioned the 
method experienced ideological developments too. Below we will discuss some of 
them. 



10.1 Light-cone sum rules 

This formalism was designed to overcome difficulties in the traditional sum rules 
for three-point functions. Assume we want to calculate a cubic coupling constant 
corresponding to the amplitude A ^ B + C, where A, B and C stand for either 
hadrons or external currents, say, electromagnetic. Typical examples are D* — > Dtt 
or 7* + 7* — i> TT (here 7* denotes a virtual photon). In many instances in order 
to get the desired constant it is necessary to apply borelization with respect to two 
momenta, pa and ps, independently. Then the components of the third momentum, 
pc-i can not be small compared to p\ b- Indeed, Pc{pa+Pb) =Pa — 'P% ^ Pa b- The 
standard condensate expansion contains a series of operators with derivatives which, 
eventually, give rise to the expansion parameter of the type PciPA + Pb) Ipa b ~ 1- 
In other words, all terms in this subseries must be summed over. 

This (partial) summation is carried out automatically if, instead of the three- 
point function {jA,jB,jc) one considers the correlation function of the currents Ja 

^"^The l/niQ expansion of the SVZ sum rules per se dates back to the work of Shuryak | |75[ ]. 
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and Jb sandwiched between the vacuum and the state \C). The vacuum expecta- 
tion values of local operators (condensates) in the SVZ sum rules are substituted by 
the light-cone wave functions {distribution amplitudes in the American literature) 
describing the momentum fraction distribution of the components in the Fock wave 
function of the hadron C. If the SVZ approach is based on the short- distance expan- 
sion of the current T products in terms of local operators, and the expansion runs 
in dimensions of the local operators, the light-cone sum rules exploit the expansion 
in nonlocal "string" operators on the light-cone; the expansion runs in twists, not 
dimensions. From the viewpoint of the standard SVZ sum rules, in the light-cone 
approach one performs a partial summation of an infinite chain of operators of ar- 
bitrary dimension, but given twist. In effect, the light cone sum rules combine the 
SVZ approach with the technique used in the description of hard exclusive processes 
(see |^8|]). The price we have to pay for the partial summation is rather high: the 
large distance dynamics is parametrized not by numbers, as in the SVZ condensates, 
but by functions - the leading twist, next-to-leading twist, and so on distribution 
functions. 

A Russian proverb says that it is better to see once than to hear hundred times. 
Let us have a closer look at a specific example. Below I will outline the calculation of 
the D*D7r constant. This example, as well as the explanations below, are borrowed 
fromRef. [^. 

The D*D7i coupling constant g is defined as follows: 

{D*{p)7fiq)\Dip + q)) = -gq,e^, 

where the hadrons' momenta are indicated in the parentheses, and is the D* 
polarization vector. To determine this constant from the light cone sum rules one 
considers the correlation function 

F^.iP,q) = i J d^xe'P-{7i{q)\T{d{xMx), c(0)z75w(0)}|0) . (84) 

The first operator in the T product is the interpolating current for D* and the 
second is the interpolating current for D. 

The correlation function (|8^ depends on two variables, and (p + qY (note 
that g2 = since the pion is on mass shell). The contribution proportional to the 
D*Dtt amplitude has two poles, in and {p + g)^, corresponding to the ground 
state mesons both in the vector and pseudoscalar channels. Namely, 

Fipq) = M!,Mn^Mn* 

m^{p^-Ml,){{p + qY-Ml) + ^"'^ 

The ellipses on the right-hand side denote a term proportional to p^ irrelevant for 
our analysis, as well as terms with single pole and non-pole continuum contributions, 
that will be exponentially suppressed under double borelization (with respect to p^ 
and {p + qY). All other notations are self-explanatory. 
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Figure 20: Light cone sum rule for 0*0^ in the leading order. The current labeled 
by p + g interpolates D while that labeled by p interpolates D*. 



On the theoretical side we assume that — M"^, and (p+g)^ — M|, are Euclidean 
and large enough to allow for the operator product expansion of the correlation 
function (P^). As usual, the corresponding Borel parameters will lie in the window. 
The simplest graph pertinent to OPE we have to build is depicted in Fig. 20. The 
large distance part of the process is denoted by the shaded area. The light quark 
lines form the pion, rather than the quark condensate, that was discussed in Sect. 
4. More exactly, we deal here with the leading twist pion wave function whose 
definition will be given shortly. 

If ^ is the fraction of the pion momentum carried by the d quark, the graph of 
Fig. 20 gives rise to the following expression: 

F, ^ q, F{p\ {p + qf) = mj^q, £ ^[pl + - ' (§6) 

where 0,r(O is the distribution function, and the ellipses denote corrections, both 
perturbative (the expansion) and non-perturbative (due to higher twists). The 
c quark Green's function is substituted, in the leading approximation, by that of a 
free quark. Emission of gluons from this line would result either in corrections 
or in the distribution functions of non-leading twists. 

For very small values of q one could have substituted the denominator in Eq. (|86|) 
by ml — p^. Then the integral / d^(f)j,{^) — >• 1, and would reduce to the matrix 
element of the operator d'j^'y^u sandwiched between the pion state and vacuum. 

Expanding the denominator in Eq. ( pUf ) in q we readily identify higher dimension 

operators that are summed over: the term linear in q originates from d Da 7/^75^, 
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quadratic in q from d DaDp Jfi^JbU, etc. The actual expansion parameter is 

pq 



p2 _ ^2 



1. 



The (leading twist) light-cone wave function is formally defined as 



{7r{q)\d{xh,^,exp{tg A„(t/)rfi/jM(0)|0) = -^qj^l d^e'^'^^MO ■ (87) 

It parametrizes the large-distance dynamics of the du pair. 

If we do the double Borel transformation, with respect to p^ and {p + qY the 
denominator in Eq. (|86D becomes 



exp(-mg .^2.^2 IM^- 



where Mi 2 are the corresponding Borel parameters, that are expected to lie close to 
each other inside the window. The sum rule is obtained by matching the double Borel 
transform of Eq. ( pSj) with the double Borel transform of Eq. (p6|). Thus, we see 
that the theoretical part of the sum rule is determined by 0,r(O.5). Of course, there 
are corrections of higher order in as and from higher twist distribution functions, 
but as usual inside the window they should be kept at a modest level. In addition 
we must check the stability of the result inside the window. Provided all this is 
properly done, one obtains |^ g in terms of (/)7r(0.5). The latter quantity must be 



extracted from independent sources. 

(The situation we face here is not quite generic. A priori one could expect the 
coupling g to be expressible in terms of (p-n^i), rather than one constant, 0^(0.5).) 

Information on 07r(^) is quite abundant; there is some limited information on the 
higher- twist wave functions too. The issue of 07r(O is almost as old as the SVZ sum 
rules themselves. The key element of the corresponding theory, which was studied 
in great detail in connection with hard exclusive processes, is the (approximate) 
conformal symmetry of QCD. Other elements are provided by the sum rule and 
lattice calculations and models. The issue is not completely settled, although it 
seems that a heated debate of the 1980 's gradually wanes. Out of two competing 
scenarios - a double hump distribution function of Chernyak and Zhitnitsky 



and a narrow distribution of Radyushkin and collaborators pO[ that is close to the 
asymptotic form of 07r(O ^ the second is gaining more recognition now. Experts 
lean towards a narrow almost asymptotic distribution. As usual, the best strategy 
is sacrificing one of the sum rules in order to extract (p-niO from data. A recent 
attempt based on the 7* -|- 7 — >■ 7r° data has been undertaken in Ref. 0]. I will 



not go into further details here, referring the reader to a very rich original literature 
devoted to this topic. 

It remains to be added that (i) the term "light-cone sum rules" first appears in 
(ii) the idea of the approach dates back to Refs. |Q ; and (iii) a review of this 



topic, with a representative list of references was published recently |8J . 
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If the partial summation of the higher-dimension operators turns out to be so 



successful in the light-cone sum rules, it is natural to ask why it is not applied 
in the vacuum correlation functions appearing in the SVZ sum rules. In this case, 
in order to partially sum, say, the quark operators we must substitute the quark 
condensate {qq) by a " string expectation value" in the vacuum, 

{q{x) exp{z^ r AMdya}q{0)) , (89) 





which is a function of x. 

In the light cone sum rules one can ascribe certain twist to each given distribu- 
tion function. Each distribution function sums up a tower of the local operators, 
with all dimensions but given twist. Therefore, the distribution functions can be 
systematically ordered in twist. Twist is obviously irrelevant in the vacuum cor- 
relation functions. Therefore, the partial summation implemented by the "string 
expectation values" (|89|) is hard to justify theoretically. A principle that would 
allow us to order distinct "string expectation values" has to be elaborated. 

Another major difference with the light cone sume rules is that the pion distri- 
bution function is known at least at a certain level, while next to nothing is known 
about the vacuum function (^91). Until recently even the large x asymptotic behav- 
ior of this function was erroneously assumed to be Gaussian. We are at the initial 
stage, when first observations are being made regarding properties of the function 
(^) following from the general structure of QCD [Q. I do not rule out that, as our 
understanding progresses, the corresponding modification of the sum rule approach 
will be elaborated. 



10.2 Heavy flavor sum rules 

Although the character of problems in which this formalism is instrumental is some- 
what different from the classical applications of the SVZ sum rules, the basic tech- 
niques - OPE, matching of the theoretical expansions with the phenomenological 
expressions represented by dispersion integrals over the observable spectral densi- 
ties, etc. - are the same. The vacuum condensates are replaced by the expectation 
values of local operators over the heavy flavor states, say, B mesons. All questions 
that can be asked in connection with the vacuum condensates exist here too; the 
answers are more transparent, however. We will consider below a sample applica- 
tion. Our main goal is discussing subtle nuances associated with the condensate 
expansion (Wilsonean OPE versus practical version, see Sect. 6). 

The heavy flavor sum rules were engineered to treat the inclusive heavy 



flavor decays, for instance the semileptonic decays B X^iv where X^. is arbitrary 
hadronic state containing c quark. The lowest-lying state is the D meson, then 
comes D* and then higher excitations/continuum. By varying the momentum of 
the lepton pair one can measure various spectra in the transition 6 — >■ c induced by 
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Figure 21: The operator product expansion for hab- For large negative e the c quark 
is far off shell, and the points of the current emission and absorption are close. The 
B meson in the initial and final states is assumed to be at rest. 

the weak current. The object of the theoretical study is the transition operator 

fabiq) =ijd'x e"^^ r{jt(x)j,(0)} , (90) 

that will be eventually sandwiched between the states \B) and {B\ (see Fig. 21). 
Here ja denotes a current of the type cVJ) with an arbitrary Dirac matrix F^; q is 
the momentum carried away by the lepton pair. 

The average of T^^ over the heavy hadron state B with the momentum pb rep- 
resents a forward scattering amplitude (the so-called hadronic tensor), 

K,{PB.q)^^^{B\fab\B) . (91) 

The structure functions Wab are obtained from the hadronic tensor by taking its 
imaginary part, 

Wab = (1/0 disc hab ■ 

All observable spectra are represented as certain integrals over Wab- The hadronic 
tensor hab and the structure functions Wab can be expanded in terms of various 

kinematic structures. For each set of the currents generically we may have up to five 
distinct kinematic structures; these details need not bother us here. For illustrational 
purposes it is sufficient to limit ourselves to the axial currents, so that Va,b lulbi 
and to the first function, 

h^v = -hig^jy ... , wi = (1/i) disc hi . 

The lowest- lying state produced by the axial current from B is D*. To end up with 
the kinematical aspects we must specify the lepton pair momentum q. It will be 
assumed that the spatial momentum g is small but non- vanishing, A <^ |g| ^ Md- 
The terms 0{(f') will be kept while those of higher order in |g| will be neglected. 
Since g is fixed, /if^ becomes a function of one (complex) variable, q^. Instead of qo 
it is more convenient to work with e defined as 

e = qomax - qo (92) 
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where 

qOmax = MB-ED*, ED'=MD*+^rf^. (93) 

ZlVl D* 

When e is real and positive we are on the physical cut where the imaginary part of 
h^^ is measurable. For negative e we are below the cut, in the Euclidean domain, 
where the T product (^) can be computed as an expansion in l/m^fe. The leading 
operator in this expansion is bb. It has dimension three. No operators of dimension 
four exist. At the next level, of dimension five, there are two operators, b{iD)'^b and 
{ig /2)lrj^^'^G fj^pb. The first operator (sometimes called the kinetic energy operator) 
is not a Lorentz scalar. In the problem considered the operators to be retained in 
the expansion need not necessarily be Lorentz scalars since, unlike the vacuum case, 
the very presence of the B meson in the initial/final state singles out a reference 
frame. 

The whole procedure is perfectly analogous to the expansion of the polar- 
ization operator discussed in Sect. 4. The operator bb occupies the same position in 
the hierarchy as the unit operator in the expansion of the vacuum correlators. The 
operators of dimension five and higher generate l/rric^b corrections analogous to the 
condensate corrections. In the problem at hand all local operators appearing in the 
expansion of the T product ( pO]) are averaged over the B meson state. In this way 
we get a theoretical expression for hf^^ off the cut. 

A bridge between h"^^ at negative e and at positive e is provided by the 
dispersion relation. Expanding the denominator in the dispersion relation in 1/e we 



obtain an infinite set of the sum rules. The third of them is as follows p6| : 



^|dee'wt%) = l^^lv' + .... (94) 

Here v is the velocity of the produced heavy hadron (in the B rest frame), and /x^ 
is defined as the expectation value 

= ^{BmD)%\B) . (95) 

This parameter, fi^, has the meaning of the average spatial momentum squared 
of the heavy quark inside B meson, and is one of the fundamental parameters of 
the heavy quark theory. One encounters fi"^ in a large number of applied problems. 
It is fair to say that in the heavy quark physics fi^ plays the same role as the gluon 
condensate in classical applications of the SVZ method. 

The ellipses in Eq. (|94D denote corrections suppressed by powers of l/m^b. 
Since the weight function in Eq. is proportional to e^, the "elastic" B ^ D* 
transition drops out, and the integral on the left-hand side is saturated by "inelastic" 
transitions, i.e. decays of the B meson in the excited Xc states. For such transitions 
Wi^{e) is proportional to ■u^ at |w| <C 1. It is convenient then to introduce a spectral 
function a{e) which does not vanish in the limit \v\ —>■ 0, 

a(e) = |t7rV^(e), e>0. (96) 
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Figure 22: A sketch of e^(T(e) in the b ^ c transition at small \v\ (solid curve). The 
spectral density a{e) is defined in Eq. (0). The dashed curve represents e^cT(e) in 
perturbation theory at one loop. 

Our brief excursion in the heavy fiavor sum rules is almost over. In terms of the 
new spectral density the sum rule we are interested in takes the form 



As with any sum rule it can be read in both directions: if /i^ is known this is a 
prediction for the spectral density. On the other hand, if the spectral density is 
measured we can take the integral and get /x^. Below we will pursue the latter 
strategy assuming that a{e) is measured. In actuality this is not the case, but we 
have a pretty good idea of how cr(e) will look like when it is measured. First of 
all, it is positive-definite. At small values of e, of order of a few units time A, 
there is one, perhaps, two resonances; then at higher values of e the spectral density 
must approach its perturbative asymptotics which is known too: cr(e) Cag/e, 
where C is a constant which was calculated in Ref. |T^. The numerical value of 
C is unimportant for our purposes. The asymptotic regime is achieved through 
oscillations - the reasons are the same as in Sects. 2.2 and 4. A schematic plot of 
e^cr(e) is depicted in Fig. 22. 

If you think about Eq. ( pTf ) you will, probably, find it remarkable: it gives a 
condensate-like non-perturbative parameter in terms of an integral over the observ- 
able spectral density. Further examination raises a question. Indeed, if the integral 
extends to infinity, it badly diverges. 

One need not be concerned, however. This infinity refiects the fact that the 
proper definition of requires a normalization point, just in the same way as the 
for the gluon condensate. In the case at hand this normalization point has a very 
transparent physical meaning - is the maximal excitation energy of the hadronic 




(97) 
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state belonging to the class of "soft excitations" (Fig. 22). Higher excitations 
are considered to be hard, and are well approximated by perturbation theory. I 
remind that the hard contributions do not belong to /i^; rather, they determine 
the coefficient functions in OPE. In the perturbative calculation of the coefficient 
functions one must remove the part with e < fi. In fact, the cut off in the excitation 
energy as a normalization point was suggested in the given context long ago in Ref. 



| 87| . Thus, in the Wilsonean approach /^^(yu) is determined as 

^^'M = ^ I' dee' a{e) . (98) 

In practical OPE we deal, instead, with subtracted "surrogates". Say, if one 
limits oneself to one-loop accuracy, 

3 



3 /"^ 

/^^" = TT / ^^^^ [^(^) " ^i-ioop(e)] 
/vr JO 



The right-hand side in this expression is fi independent at one loop. Its fi dependence 
shows up only at two loops. If we agree to discard all effects 0{al) and higher, 
just totally ignore them, then the integration in "/z^" can be extended to infinity. 
Correspondingly, the coefficient functions must be calculated at one-loop, following 
the standard Feynman rules, with no removal of the e < /x part. 



11 Sum Rules and Lattices 

In this Lecture it has been already noted, more than once, that certain aspects 
of the condensate approach, as well as a rich experience in the analysis of various 
spectral densities, might be useful for lattice practitioners. Remarks to this effect 
are scattered here and there. Of particular importance are qualitative observations, 
such as an enhanced sensitivity to the vacuum fermion loops in certain channels (this 
effect is neglected in the quenched approximation), the existence of new, abnormally 
large scales in the 0^ glueball channels, and so on. Surprising though it is, all these 
hints, which might have produced a strong impact on the lattice theory, are so far 
almost totally ignored The reaction of rejection of any ideas coming from outside 
the lattice theory per se plagues the lattice community. This happens even in those 
cases where the overlap is quite obvious. 

Let me give an almost anecdotal example. As was mentioned in Sect. 7, in 
1981 a wide class of "QCD scale anomaly" low-energy theorems was obtained [Q] in 
connection with the SVZ sum rules. These theorems relate to each other n-point 
functions with arbitrary number of insertions of the operator (y{x) = ^^(x) oc G' 
(see Eq. (^9])) at vanishing momentum. One of this theorems can be cast in the 

^^I am aware of two dedicated works where the sum rule resuhs where analyzed in parallel with 
the lattice calculations [p8| . This seems to be a rare exception. 
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Figure 23: Rectangular contour in the Wilson loop operator. One can consider it 
as a trajectory of infinitely heavy quarks (in the imaginary time). R is the distance 
between the quark and antiquark, T — > oo. 



form of a prediction for the expectation value of in the presence of the Wilson 
loop operator W, 

{W)-' Q d''xa{x)W^^ = V{R) + R^^^ (99) 

where V{R) is the static potential between the heavy quarks separated by distance 
R. 



W = <j> Tr ex^{ig A^{x)dx 



the contour C in the definition of the Wilson loop is depicted on Fig. 23, and the 
subscript E reminds us that all definitions and calculations refer to the Euclidean 
space. Moreover, the integration in Eq. (^) is performed in three directions per- 
pendicular to the Euclidean time; the subscript c refers to the connected part, i.e. 
the disconnected part of the correlation function (|99D must be discarded. One can 
interprete the operator a{x) in Eq. ( P5| ) as a probe of the energy density of the vac- 
uum medium at the point x in the presence of a pair of heavy static quark sources 
separated by distance -R 0. In the case of purely linear potential the right-hand 
side reduces to 2V{R). In this form the theorem (|99D was derived in Ref. (see 
Appendix B) . The generalization to the case of arbitrary potentials is due to Dosch 
et al. |8|. 

Derivation of the theorem (|99|) is very similar to other theorems discussed in Sect. 
7. The functional-integral representation for (W) is differentiated with respect to 
I/Qq. We then use the fact that 

The latter equality, in turn, follows from consideration similar to that presented in 
Sect. 7. Thus, the low-energy theorem ( p9D is a reflection of the scale anomaly of 
QCD, and as such, represents a fundamental aspect of the theory. 

Now I come to a culmination of the story. In 1987 the theorem ( ^9]) was redis- 
covered in the lattice community |^0|) where it goes under the name of the "action 
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Michael sum rule". The original derivation was erroneous, the second term on the 



right-hand side of Eq. (pQ]) was omitted. The error was corrected in Ref. |89 



which was specifically devoted to the issue. Nevertheless, the theorem continues to 



circulate in the lattice community as the "Michael sum rule" , see e.g. pT |. 

This example of a "lattice xenophobia" is by no means unique. Needless to say 
that I consider the situation as absolutely unhealthy. The problems we face in the 
hadronic world are too complicated to afford neglecting the work done in related 
areas. The analytic and numerical approaches should complement each other. I 
sincerely hope that the attitude of the lattice community will start changing after 
recent breakthrough discoveries in supersymmetric gauge theories. Supersymmetry 
is a powerful tool which allows us to reveal features of strong coupling gauge dynam- 
ics we never suspected of. One of the lessons which seems undeniable is the proof 
that massless (or light) quarks play much more significant a role than is universally 
believed among the lattice practitioners. In the quenched approximation, the most 
wide spread today, the dynamical quarks are neglected altogether. Even the most 
advanced lattice investigations, which do not rely on the quenched approximation, 
routinely use extrapolations in the number of the massless quarks. Supersymmetric 
gauge theories teach us that increasing the number of massless quarks from two to 
three or changing the color representation of the fermion fields from fundamental to 
adjoint is not necessarily a smooth process. Dynamics may change drastically! 

Penetration of ideas and insights obtained analytically, in the practice of lattice 
calculations, where they can be used, at the very least, to test reliability of various 
approximations, is unavoidable. First steps in this direction have been already 
reported |93 . 



12 Vacuum Fluctuations Are Subtle Creatures 

The rich SVZ phenomenology confirms the wide-spread belief that the structure of 
the hadronic family, in all its gross features and intricacies, derives from a compli- 
cated organization of the vacuum medium. The vacuum fiuctuations are violent, but 
fine-tuned in a very specific way. Within the method, as it exists today, we do not 
calculate the vacuum condensates. This does not mean, of course, that the origin 
of the condensates, and their relative values, are of no interest. I merely wanted to 
say that the problem of the vacuum structure in analytic QCD is not solved in full. 
There are indirect methods, however, which show that this structure is very subtle 
- much more subtle than one might think of a priori - and depends on seemingly 
insignificant details of the theory. For instance, changing the number of colors from 
three to two dramatically changes the pattern of the spontaneous chiral symmetry 
breaking resulting in a drastic rearrangement in the Goldstone sector of the hadron 
family. Instead of Nj — 1 QCD "pions" [Nf stands for the number of the light quark 
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flavors) we get 2Nj — Nf — 1 massless Goldstone hadrons Q. This conclusion can be 
achieved |^ from consideration of the 't Hooft matching conditions (for earlier 



analyses based on less rigorous arguments see Ref. 

Let us dwell on this issue. Assume that we have Nf massless (light) quarks 
but the color gauge group is SU{2). The quarks belong to the fundamental rep- 
resentation of SU{2). Unlike QCD, in which the color triplets and antitriplets are 
independent representations, the SU{2) doublets are essentially the same as anti- 
doublets. In other words, all representations of SU{2) are (quasi)real. 

This simple observation has far reaching consequences. Indeed, the flavor sym- 
metry of the Lagrangian is SU{2Nf) now, rather than SU{Nf) x SU{Nf) we deal 
with in QCD. Out of 4Nj — 1 currents 2Nj — Nf — 1 are axial; all correspond to 
the spontaneously broken symmetries since there is no way one can achieve the 
matching of the triangular AVV anomalies otherwise ||9^. This means that the 
pattern of the chiral symmetry breaking is SU{2Nf) — > Sp{2Nf) rather than the 
SU{Nf) X SU{Nf) SU{Nf) we got used to in QCD. The very gross features of 
the vacuum structure do depend on the number of colors! 

Not only do they depend on the number of colors, they are also sensitive to the 
fermion contents of the theory. 

Indeed, let us substitute the standard quarks of QCD with one Majorana fermion 
A in the adjoint representation of the color group. The theory thus obtained is 
nothing but supersymmetric gluodynamics. Superficially it looks very similar to 
conventional QCD. At first sight, one can hardly expect any conspicuous deviations 
from the conventional pattern of behavior. 

It has been known for many years that the supersymmetric gauge theories posses 
some miraculous properties. In particular, in supersjTiimetric gluodynamics the 
Gell-Mann-Low function is known exactly [RBI 



- 1 - (Na./2n) ' f^"^' 

where the SU{N) gauge group is assumed. 

As QCD, supersymmetric gluodynamics is believed to be confining in the in- 
frared domain. Unlike QCD, however, the condensate (A^) may or may not develop. 
It was argued recently [0 that two distinct phases coexist in this theory - the con- 
ventional chirally asymmetric phase with (A^) 7^ 0, that reminds QCD, and a very 
unusual chirally symmetric phase where (A^) = 0. The arguments leading to this 
conclusion are too technical to be reproduced here. They would lead us far astray. 
The interested reader is referred to Ref. p7| . 

Instead, let us dwell on another example which is pretty close to actual QCD. In 
fact, we will consider just conventional QCD with more than three massless flavors. 



^^It is curious to note that the existence of the extra massless states was discovered in the lattice 
SU{2) theory, which is being used as a theoretical laboratory for decades, only recently. The result 
is known in analytic QCD since the mid-eighties. This is another manifestation of the "lattice 
xenophobia", see Sect. 11. 
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The Gell-Mann-Low function in QCD has the form 



2n 47r2 

2 19 
b = 11 --Nf, bi = 5l-—Nf. (102) 

At small as it is negative (asymptotic freedom!) since the first term always domi- 
nates. With the scale fi decreasing the running gauge coupling constant grows, and 
the second term becomes important. Generically the second term takes over the first 
one at a^/vr ~ 1, when all terms in the expansion are equally important, i.e. in 
the strong coupling regime. We do not know what exactly happens at as/vr ~ 1. As- 
sume, however, that for some reasons the first coefficient b is abnormally small, and 
this smallness does not propagate to higher orders. Then the second term catches up 
with the first one when a^/vr -C 1, we are in the weak coupling regime, and higher 
order terms are inessential. Inspection of Eq. (|102|) shows that this happens when 



Nf is close to 33/2, say 16 or 15 {Nj has to be less than 33/2 to ensure asymptotic 
freedom). For these values of Nf the second coefficient bi turns out to be negative. 
This means that the /3 function develops a zero in the weak coupling regime, at 

^ = -^«l. (103) 

ZTT —bi 

(Say, if Nf = 15 the critical value is at 1/44 and is indeed small.) This zero is 
nothing but the infrared fixed point of the theory. At large distances — > a*, and 
/9(q;*) = 0, implying that the trace of the energy-momentum vanishes. Then the 
theory is in the conformal regime. There are no localized particle-like states in the 
spectrum of this theory; rather we deal with massless unconfined interacting quarks 
and gluons; all correlation functions at large distances exhibit a power-like behavior. 
In particular, the potential between two heavy static quarks at large distances R 
will behave as ~ <y*/R, i.e. we have a pure Coulomb behavior. The situation is not 
drastically different from conventional QED. As long as a* is small, the interaction of 
the massless quarks and gluons in the theory is weak at all distances, short and large, 
and is amenable to the standard perturbative treatment (renormalization group, 
etc.). The chiral symmetry is not broken spontaneously, the quark condensates do 
not develop. QCD becomes a fully calculable theory. 

The fact that at Nf close to 16 QCD becomes conformal and weakly coupled in 
the infrared limit is known for about 20 years P5| , . If the conformal regime takes 



place at Nf = 16 and 15, let us ask ourselves how far can we descend in Nf without 
ruining the conformal nature of the theory in the infrared? Certainly, when Nf is 
not close to 16 the gauge coupling a* is not small. The theory is strongly coupled, 
as conventional QCD, and, simultaneously, conformal in the infrared. The chiral 
symmetry presumably is unzbroken, (qq) = 0. The vacuum structure is totally 
different, and so are all properties of the theory. 
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The range of Nf where this phenomenon occurs is called the conformal window. 
The right edge of the conformal window is at Nf = 16. There are indications that 
the left edge of the conformal window A'^^ lies not too far from actual QCD (i.e. not 
too far from Nf = 3). 

First, in the instanton liquid model it was found ||100[ | that at Nf = 5 the chiral 
condensate disappears, {qq) = 0. Although nothing is said about the onset of the 
conformal regime in this work, the vanishing of (qq) may be a signal. More definite 
is the conclusion of the lattice investigation ||101|| . It is claimed that at Nf = 7 
not only the quark condensate disappears, {qq) = 0, but the decay law of the 
correlation functions at large distances changes from the exponential to power-like, 
i.e. the theory switches from the confining regime to the conformal one. One could 
speculate that at Nf = 5 and 6 we are in the intermediate phase, when confinement 
is still operative at large distances but the chiral symmetry is unbroken. Previously 
it was believed that that confinement of color in QCD implies the chiral symmetry 
breaking ||102|| . Today's wisdom tells us that this need not be necessarily the case. 
At least in some supersymmetric gauge theories we have both: confinement and 
unbroken chiral symmetry (see e.g. In any case, we see that ascending from 

Nf = 3 to Nf = 5 or 7 - quite a modest variation- we drastically change the vacuum 
structure of the theory, with the corresponding abrupt change of the picture of the 
hadron world. 

An interesting question is what happens slightly below the left edge of the con- 
formal window Nf. If Nf were a continuous parameter, and the phase transition 
in Nf were of the second order, slightly below Nf the string tension a would be 
parametrically small in its natural scale given by A^. In actuality Nf changes dis- 
cretely. Still it may well happen that at Nf = NJ — 1 the ratio cr/A^ is numerically 
small. This would mean that the string and its excitations, whose scale is set by 
a, are abnormally light. Under the circumstances one might hope to build effective 
low-energy approaches analogous to the chiral Lagrangians of actual QCD. In the 
limit (j/A"^ —>■ the string becomes, in a sense, classic; the quantum corrections are 
unimportant. I do not rule out that the string representation of QCD - the holy 
grail of two generations of theorists - is easier to construct in this limit. One could 
even dream of an expansion in NJ — 3. At this moment, this is a pure speculation, 
however, and I have to wind up. 



13 Instead of Conclusions 

As time passes the hopes that the full analytic solution of QCD will be found fade 
away. After all, the problem is with us for over a quarter of the century, and none 
of numerous theoretical attacks reached the goal. Many theorists whose philosophy 
is "all or nothing" abandoned the field Does it mean that the theory of hadrons 

^^The "all or nothing" philosophy is wide-spread and, unfortunately, not only in theoretical 
physics. This is a favorite child of the so called revolutionaries in all times and in all countries. 
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is an aging science at the verge of submerging into a permanently dormant state? 

I do not think so. Those theorists who stayed in the field can claim many partial 
successes. The field is messy and down-to-earth but - what can we do? - this is the 
only world God gave to us. The theory of hadrons is alive. Our understanding of 
the vacuum structure of quantum chromodynamics and various dynamical regimes 
it can support continues to grow. Especially fruitful were the last years, with many 
breakthrough discoveries in supersymmetric gauge theories. These discoveries are 
to be used as hints and insights in actual QCD. 

OPE-based methods - and the SVZ sum rules is one of them - continue to grow 
too. They are admittedly approximate, but their analytic nature and a graphic phys- 
ical interpretation make them indispensable in many practically important problems. 
In some issues they share their role with other methods, which were developed later, 
for instance, lattice QCD. Since the theory of hadrons is so difficult it is in our best 
interests to combine all sources of information. We are in no position to neglect one 
of them in favor of the other. Many ideas and theoretical devices surfaced in con- 
nection with the SVZ method (e.g. the low-energy theorems. Sect. 7); the method 
also produced a number of variations and spin-offs. This is a healthy process which 
will hopefully continue. 
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The misfortunes it brought to our world are innumerable. Needless to say, I personally think this 
is a rotten philosophy. 
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14 Appendix. Some Useful Definitions and For- 
mulae. 



The Gell-Mann-Low function is defined as 

das 



(91n/i 



27r 47r2 



+ ... 



(A.l) 



This definition is standard except that in some sources the first coefficient is called 
Po, the second and so on (cf. Ref. P). At two loops the /3 function is scheme 
independent. 

The coefficients are 



2 19 
n--Nf, bi = 51-—N 



f ■ 



(A.2) 



The running coupling as{fi) is parametrized as follows 



An I ^ 2bi In In 



Mn^ 



(A.3) 



A more convenient form of the very same expression is 

27r , ^1 1 1 /^^ , 

— -— = In — + — Inln — + ... 

asifi) A b A2 



(A.4) 



Being expanded, at two loops, the latter expression is equivalent to the former; Eq. 
( |A.4| ) effectively sums up some higher order terms. 
The scale parameter A following from Eq. ( [A. 3D is 



A^ 



2\h/b 



/i''exp 



27r 



27r 



-^In 



(A.5) 



The first factor on the right-hand side is an (inconvenient) artifact of the definition. 
It could have been easily avoided. Since the definition is standard, we will keep it 
not to cause confusion. 

Perturbative calculations in QCD, as a rule, are carried out in the so called 
modified minimal subtraction (MS) scheme ||103|| , using dimensional regularization. 
For massless quarks and gluons it works nicely; the treatment of the heavy quark 
mass thresholds in this approach is ugly. The most common is the step-function ap- 
proximation ||104|| : immediately above the threshold the quark is declared massless, 
while below the threshold it is treated as infinitely heavy, so that it is frozen out and 
does not participate in loops. This is equivalent to treating Nj as a function of fi 
constructed from several step functions, Nf = 3 below rric, then it jumps to Nf = 4 
between rric and nib, etc. Matching conditions at thresholds require equivalence of 
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one effective theory with Nf massless quarlcs to another effective theory with Nf — 1 
massless quarlcs. 

If the step-function approximation is apphed, at one loop the couphng constant 
is obviously continuous, while the first fi derivative experiences a jump. At two and 
three loops the coupling itself becomes discontinuous if the matching is done at the 
quark masses ||105|| . Several suggestions as to how one can smooth out the running 
coupling constant using various "physical" definitions were presented in the litera- 




ture 



In 



I propose a somewhat different approach inspired by supersymmetry. 
supersymmetric theories the quark mass thresholds can be accounted for exactly, to 
all orders ||107|| . For instance, consider supersymmetric QCD with one flavor, with 
the mass term m. Assume that we start our evolution at a high normalization point 
Mq (the corresponding coupling constant is a^o), pass the threshold and descend 
down to /X <^ m. A continuously running ag can be obtained from the formula 



2tt _ 2n Mo ^ Mq 



(A.6) 



for SU(N) gauge group. Here Z{fj,) is the Z factor of the matter fields and fj, is 
arbitrary: larger or smaller than m. At n = m and below the second logarithm 
freezes at InMo/mo where mo = mZ is the mass parameter normalized at Mq. The 
one loop Z factor implies a;s(/i) at two loops. Needless to say that Z{fi) varies 
continuously at one loop. At the two-loop level the high- and low-energy scale 
parameters are related as follows: 
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where 
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(A.7) 



2N 



At this level of accuracy mo and a^o in the square brackets can be treated in the 
leading logarithmic approximation. In this approximation the expression in the 
square brackets is renormalizat ion-group invariant. In principle, the procedure can 
be extended to all loops, but we will not pursue this goal here. 

In non-supersymmetric QCD the exact treatment of the mass thresholds seems 
impossible. However, numerically the situation is quite close to what we have in 
supersymmetric QCD. One can work out compact formulae applicable at two loops. 
Say, for the charm threshold 
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and 



Correspondingly, 
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If it were not for the last factor (which is quite close to unity) the relation between 
the low- and high-energy scales would be similar to that in supersymmetric QCD. 
Matching at two thresholds, charm and beauty, yields 
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